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Abstract 

I describe the history of the calculations of NLO and NNLO QCD corrections to weak decays 
of mesons and particle-antiparticle mixing in the period 1988-2010. Also existing calculations 
of electroweak and QED corrections to these processes are included in this presentation. These 
efforts bear some analogies to the climbing of Himalayas and various expeditions by several teams 
of strongly motivated "climbers", acting mainly in Germany, Italy, Poland and Switzerland, 
allowed to move this field from the LO through the NLO to the NNLO level. The material 
is meant to be a guide to the rich literature on NLO and NNLO corrections in question and 
includes several anecdotes related to the climbs that I was involved in. While, by no means this 
presentation should be regarded as a review of the size of NLO and NNLO corrections to weak 
decays, I hope that some of the comments made in the course of the presentation could turn out 
to be not only amusing but also instructive. 



* Dedicated to the members of the Munich NLO Club in Table 1. 
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1 Introduction 



In April 1988 the workshop on "Hadronic Matrix Elements and Weak Decays" took place in the 
Ringberg castle at the Tegernsee lake near Munich. This workshop, organized by Jean-Marc 
Gerard and myself, with a great help of Willy Huber, the last secretary of Werner Heisenberg, 
gathered a large portion of experts that were making the first steps towards the calculation 
of hadronic matrix elements relevant for — mixing, mixing and non-leptonic 

decays of K mesons, in particular K — t- vrvr. Also non-perturbative aspects of semi-leptonic K- 
meson decays belonged to the important topics of this workshop. The representatives of lattice 
calculations, large N approach, QCD sum rules, hadronic sum rules and chiral perturbation 
theory were presenting their views on the subject. In particular Bill Bardeen summarized 
the large N approach to weak non-leptonic K meson decays developed by him, Jean-Marc 
Gerard and myself in 1986-1988 [H [21 [3l H]. Using this approach we were able to obtain the 
first quantitive, even if approximate, results in QCD for the matrix elements relevant for the 
A / = 1/2 rule, e' /e and — mixing (Bk). These results were certainly not appreciated 
by other groups during this workshop that were confident to obtain much better results in the 
following years. In particular Luciano Maiani and Guido Martinelli promised us to provide in 
a few years a much better explanation of the observed A / = 1/2 rule in K ^ vrvr decays 
within the lattice framework as well as the value of I3k in QCD. Also representatives of chiral 
perturbation theory and hadronic sum rules were rather critical about our work. But in 1988 
only very few understood our approach, the whole field was in its infancy and it is not surprising 
that all competing groups had rather differing views on the subject. 

The workshop was certainly a great success with hot discussions in essentially all rooms of 
the castle and several participants not leaving it for the five days of the workshop. In spite of 
this it was rather clear to me that I do not want to take part in this enterprise any longer. I was 
sceptical that one could improve the accuracy of our calculations and I certainly did not want 
to be a member of a big lattice group or joining QCD sum rule and chiral perturbation groups 
that were active already for many years. 

During the last supper of the Ringberg workshop Guido Martinelli and me realized that it 
would be important to calculate NLO QCD corrections to the Wilson coefficients of penguin 
operators relevant for K — )• vrvr decays. In 1981 Guido took part in the pioneering calculation of 
the two loop anomalous dimensions of the current-current operators. This calculation done in 
collaboration with Guido Altarelli, Curci and Petrarca O |6] has been unfortunately performed 
in the dimensional reduction scheme (DRED) [7J that was not familiar to most phenomenologists 
and its complicated structure discussed in detail by these authors most probably scared many 
from checking their results. Moreover it was known that the treatment of 75 in the DRED 
scheme, similarly to the dimensional regularization scheme with anticommunicating 75 (known 
presently as the NDR scheme), may lead to mathematically inconsistent results. Consequently 
it was not clear in 1988 whether the result of Altarelli et al. was really correct. Being a member 
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of the theory group of the Max-Planck Institute for Physics in Munich (MPI) for six years I was 
exposed very much to this problematic. Peter Breitenlohner and Dieter Maison [8] were rather 
critical about the DRED and NDR schemes. According to them and other field theorists only 
the 't Hooft-Veltman (HV) scheme [9] for 75 was mathematically self-consistent. However, this 
scheme was also not familiar before 1988 to many phenomenologists. 

While visiting Technical University of Athens in 1984 I learned about the second two-loop 
calculation of anomalous dimensions of current-current operators. Two young greek physicists, 
Tracas and Vlachos [10], performed in 1982 the Altarelli et al. calculation of 1981 in the NDR 
scheme, obtaining the result that differed considerably from the one of the Italian group. They 
could not clarify the reason for this discrepancy and in 1984 being involved heavily in other 
projects I simply did not have time to have a closer look at this problem. 

At this last supper of the Ringberg 1988 workshop Guido told me that he will put some 
of his PhD students to look into NLO QCD corrections to Wilson coefficients of QCD penguin 
operators relevant for K — )• vrvr decays and I told him that I will look at this problem as well. 
However, in April 1988 I was still at the MPI and did not have any PhD students who could join 
me in this enterprise. Moreover, due to heavy involvement in the organization of the ICHEP 
1988 in Munich and other time consuming matters like the proceedings of the 1988 Ringberg 
workshop [11], lectures on our large N approach to weak decays at the Summer School in Jaca 
(Spain) |12j and most importantly because of my moving from MPI to the Technical University 
Munich, I did not have time to start this new project until October of the same year. 

During my summer vacation 1988 I read several books about the Himalaya expeditions. 
Among them the one by colonel Hunt, in which he described in detail the well known 1953 
Mount Everest expedition that he organized. From these books I learned also about the compe- 
tition between Reinhold Messner and polish climber Janusz Kukuczka to conquer the fourteen 
highest Himalayan summits, the ones over 8000 m 0. These were truly fantastic achievements 
but I wondered whether the difficulty of climbing a 8000 m high mountain by an experienced 
mountainer could be comparable in 1988 to the difficulty of a NLO calculation of weak decays 
performed by an experienced physicist like me. This comparison is not fully idiotic. After all 
the difficult pioneering NLO and NNLO calculations in the last 22 years required not only high 
technical skills but also certain planing in advance and first of all psychological and physical 
strength. I mean here the ability to be involved in a calculation that results sometimes in a 
single number but lasts at least six months and often a year or longer. The air during NLO and 
NNLO calculations can be very thin indeed. 

These thoughts prompted me to generalize my plan for the NLO analysis of — t- vrvr decays 
to all relevant K and B decays including rare, radiative and in particular CP-violating decays. 
In 1988 this field was, with respect to NLO QCD corrections, essentially unexplored and all NLO 

^Messner won this competition but Kukuczka was the second to reach aU these highest summits by October 
1988. Unfortunately he died on October 24th, 1989 in an avalanche on the South Face of Lhotse. 
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summits were still waiting to be conquered. Being the first to complete all these calculations 
would certainly be an achievement with a lasting impact on the phenomenology of weak decays. 
These were my dreams of 1988. Feeling like colonel Hunt before the Mount Everest expedition 

I made a list of most interesting decays and the corresponding operators. This list is given in 
Section 2.2 and in various Tables below. As mentioned before only the current-current operators 
Qi and Q2 have been studied at NLO before 1988 but the status of these calculations was unclear. 

The next step was to find a team of physicists with whom I would perform these calculations. 
To do them alone in 1988 would be a pure madness. Like Guido Martinelli I could in principle 
count on PhD students but in October 1988 I had none who could be put on this project 
and even if I was hoping to get some students at TUM soon, it was not certain that it would 
happen. Moreover, the knowledge of gauge theories at TUM in 1988 was very limited (Hans 
Peter Nilles and me were hired to change this situation) and without a series of lectures on 
renormalization, renormalization group methods and loop calculations, sending young students 
to do NLO calculations in QCD would be impractical and certainly irresponsible. I estimated 
that before the fall of 1989 I could not count on any help from my future PhD students and/or 
post-docs that were supposed to arrive in October 1989. 

However, I certainly could not wait until the fall of 1989. Afterall I was convinced that 
Guido already worked on this project with his students. Therfore I told Jean-Marc Gerard, 
who was at the MPI at that time, about my plans. Between 1984 and 1988 we have written 

II papers together and I was convinced that he was the right person for this grand project. 
Unfortunately Jean-Marc did not want to join me in this expedition. He basically told me that 
I was crazy to think about calculating NLO corrections to weak decays that were polluted by 
hadronic uncertainties. Li principle I could also ask Bill Bardeen with whom I did my first 
NLO calculations for deep inelastic structure functions and photon structure functions ten years 
earlier \13\ [H] . But we were seperated by the Atlantic and moreover I had some doubts that 
Bill would be interested in this project. 

On my last day at MPI, the members of the MPI theory group were giving 5 min talks about 
their research, in the spirit of similar talks in the Theory Group at CERN. At this meeting I 
informed my MPI colleagues about my project and that I was looking for collaborators. There 
was no reaction. I left MPI rather frustrated. 

Few days later, sitting already in my office at TUM, I got a phone call from Peter Weisz, 
who joined MPI few months earlier. I knew Peter from his work with Martin Liischer, but as 
his field was rather different from mine I had only a few conversations with him until then. To 
my great surprise, Peter was very much interested in my project and asked me whether he could 
join me in this enterprise. I told him that I was delighted. On this day the Munich NLO Club 
(MNLC) was born. The club consisted only of two members but our team could start the first 
climb. 

24 years later I can report that the MNLC consists of 31 members, most of them working now 
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Table 1: Members of the MNLC. My PhD students are given in boldface and the PhD students 
of my PhD students (my grandchildren) in "slanted" . The remaining members were assistants, 
post docs or visitors in my group or at MPI. 



Patricia Ball 


Matthaus Bartsch 


Guido Bell 


Christoph Bobeth 


Stefan Bosch 


Joachim Brod 


Gerhard Buchalla 


Andrzej J . Buras 


Kostja Chetyrkin 


Andrzej Czarnecki 


Thorsten Ewerth 


Robert Fleischer 


Paolo Gambino 


Jennifer Girrbach 


Martin Gorbahn 


Ulrich Haisch 


Stefan Herrlich 


Matthias Jamin 


Sebastian Jager 


Axel Kwiatkowski 


Markus E. Lautenbacher 


Alexander Lenz 


Mikolaj Misiak 


Manfred Miinz 


Ulrich Nierste 


Gaby Ostermaier 


Volker Pilipp 


Nikolas Pott 


Emanuel Stamou 


Jorg Urban 


Peter Weisz 







outside Munich. Their names are listed in Table [T]E|. The grand project that I outlined in 1988 
and started with Peter Weisz soon after turned out to be a great success. Peter was an active 
member of the MNLC only in the period 1988 — 1992 but these were very important years and 
his participation had an invaluable impact on the full project. The project has been completed 
within the Standard Model (SM) in the first years of this millennium. Also some calculations 
beyond the SM, in particular within the MSSM, have been done. We were not always the first 
to climb a given NLO summit, but MNLC is basically the only group that calculated NLO 
corrections to all relevant decays within the SM. In the last decade several NNLO calculations 
have also been performed in our club. 

Three generations of physicists took part in this enterprise with my PhD students in Table [1] 
given in boldface and the PhD students of my PhD students in "slanted" . Bartsch, Bell, Bosch 
and Pilipp were PhD students of Gerhard Buchalla at the Ludwig-Maximilian University in 
Munich and Joachim Brod graduated in Karlsruhe under the supervision of Uli Nierste, becoming 
the member of MNLC only two years ago as a PostDoc in the group of Martin Gorbahn. Emanuel 
Stamou got his Diploma under supervision of Martin Gorbahn. The remaining members except 
for Peter Weisz and myself were assistants, post docs or visitors in my group or at MPI. 

In the days of SPIRES and INSPIRE it is easy to verify that our work had an impact on 
particle physics. Around 100 papers on NLO and NNLO corrections have been published by the 
members of the MNLC, where mainly papers are counted in which NLO and NNLO calculations 
have been performed and not papers in which only the phenomenological implications of these 

*A NLO calculation for a weak decay in Munich is necessary for a membership in the MNLC. 
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calculations have been analyzed. As you will see at the end of this writing, the papers from our 
club amount to roughly 2/3 of all papers in the field of weak decays in which such calculations 
have been performed. Our papers collected over 14 000 citations with 19 papers above 250 
citations and 23 with 100+ . Moreover the Hirsch number, taking into account that only roughly 
90 papers have been published, is rather impressive: h=57. This number will certainly increase 
in time. 

The purpose of the following presentation is the recollection of these activities and a summary 
of the present status of NLO and NNLO calculations in weak decays. I have organized the mate- 
rial in the following manner. In the next section I will summarize the theoretical framework for 
weak decays. This will be a compact presentation to which I will refer frequently in subsequent 
sections. The full exposition of the technicalities that I will try to avoid as much as possible can 
be found in the Rev. Mod. Phys. article written in collaboration with Gerhard Buchalla and 
Markus Lautenbacher [T5], my Les Houches lectures [16j and of course in the original papers. 
Section 3 is devoted to NLO QCD corrections to AS = 1 and AB = 1 non-leptonic decays. I will 
be rather detailed about the structure of QCD corrections to these decays because the operators 
involved there play an important role, even if indirectly, in essentially all weak decays. Also 
existing NNLO calculations for these decays will be summarized. Section 4 describes AS" = 2 
and AB = 2 transitions in some detail and very briefly AB = transitions. Section 5 is devoted 
to rare K and B decays, in particular K'^ — t- tt^vu, Kj^ — t- ir^uiy and Bs^d — ^ • Section 6 

is devoted to the K2 of weak decays, the inclusive decay B — )• Xsj with a few comments on the 
B ^ Xg gluon decay. In addition to a detailed description of the history of NLO calculations we 
will summarize the present status of NNLO calculations. Finally, we will list the literature for 
the corresponding exclusive decays B — )• K*{p)'y. In section 7 we discuss the NLO corrections 
to Kl 7r°/+/~ and in Section 8 the NLO and NNLO corections to B ^ Here also 

the decays B — )■ K*{p)l~^l~ will be mentioned. 

A very special Section is Section 9 because it is not written by me but by one of the prominent 
members of MNLC, Gerhard Buchalla, who is also one of the fathers of the QCD Factorization 
approach to non-leptonic decays. I thought that the QCD calculations in this approach should 
also have a place in this presentation and I asked Gerhard to help me in this matter. We conclude 
in Section 10. 

While in certain parts of our review we will enter some details, the material is meant to be 
a guide to the rich literature on NLO and NNLO corrections to weak decays. It is certainly 
not a review of the subject in the spirit of \15\ [TBj . I just wanted to summarize what has been 
done during the last twenty two years, listing in particular the first climbs of the existing NLO 
and NNLO summits and few subsequent climbs that used different methods or routes to reach 
a given summit. Thus the full material can be considered as a chronicle of NLO and NNLO 
calculations (1988-2010) in weak decays with several anecdotes behind the scene related to the 
climbs that I was involved in and several, hopefully, instructive comments for non-experts that 
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probably are hard to find in the most recent very technical literature on NNLO corrections. 

Before describing my NLO-story in more explicit terms I will make an express review of the 
theoretical framework for weak decays summarizing at the end the present status of NLO and 
NNLO calculations that are discussed in detail in the subsequent sections. I should stress that 
several NLO and NNLO QCD calculations in the framework of SCET, QCD sum rules, light- 
cone sum rules, for non-leading terms in heavy quark expansions, heavy quark effective theory, 
charmonia are left out from this presentation because I did not participate in these studies. 
With the help of SPIRES or INSPIRE clicking the names of Ball, Bauer, Beneke, Brambilla, 
Chetyrkin, Hoang, Jamin, Kiihn, Neubert, Stewart and several other masters of these fields one 
can easily find all papers on these topics. The stories behind these calculations are unknown 
to me. In this context one should mention in particular numerous papers of Matthias Jamin 
on QCD corrections relevant for QCD sum rules, numerous papers by the Karlsruhe QCD club 
lead by Hans Kiihn, in particular their results on quark masses, and the work of Andre Hoang 
and his collaborators among the others. 

The following section is rather heavy but I hope that it will facilitate the reading of the 
subsequent sections for non-experts. Experts, knowing this technology, can skip this Section to 
go directly to Sections 3-10 in order to check whether I cited them properly. On the other hand 
the classification of QCD corrections into eight classes, presented in Section 2.6 could also be of 
interest to them. 

2 Theoretical Framework for Weak Decays 
2.1 OPE and Renormalization Group 

The basis for any serious phenomenology of weak decays of hadrons is the Operator Product 
Expansion (OPE) [17\ I18j . which allows to write the effective weak Hamiltonian simply as 
follows 



Here Gi? is the Fermi constant and Qi are the relevant local operators which govern the decays 
in question. As we will see below they are built out of quark and lepton fields. The Cabibbo- 
Kobayashi-Maskawa factors V^j^j^ [191 [20] a-^d the Wilson coefficients Cj(/i) describe the strength 
with which a given operator enters the Hamiltonian. The latter coefficients can be considered as 
scale dependent "couplings" related to "vertices" Qi and can be calculated using perturbative 
methods as long as the scale /.i is not too small. 

An amplitude for a decay of a given meson M = K,B, .. into a final state F = ttvv, tttt, DK 
is then simply given by 




(2.1) 



A{M ^F) = {F\n,ff\M) = ^Y.^i.KMC^{^^){F\Q^i^^)\M) 



(2.2) 
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where {F\Qi(fj,)\M) are the matrix elements of Qi between M and F, evaluated at the renor- 
malization scale fi. 

The essential virtue of the OPE is this one. It allows to separate the problem of calculating 
the amplitude A{M — t- F) into two distinct parts: the short distance (perturbative) calculation 
of the coefficients Ci{fj,) and the long-distance (generally non-perturbative) calculation of the 
matrix elements {Qi{^)). The scale ^ separates, roughly speaking, the physics contributions into 
short distance contributions contained in Ci{iJ,) and the long distance contributions contained in 
{Qi{n)). Our presentation is mainly devoted to the calculations of the coefficients Ci{p) within 
the SM and some of its extensions. Only in the case of inclusive B decays we will also look at 
the evaluation of the actual decay as it can be done in perturbation theory. 

It should be stressed at this point that our presentation would not exist without the asymp- 
totic freedom in QCD |21l [22] that allows the calculations of Wilson coefficients by means of 
ordinary or renormalization group improved perturbation theory. The precision of these calcu- 
lations increased in the last twenty years not only because of NLO and NNLO QCD calculations 
but also because of the more accurate determination of the strong coupling as for which the 
most recent result reads [23]: 

af^{Mz) = 0.1184 ± 0.0007 . (2.3) 

Now, the coefficients Cj include in addition to tree level contributions from the VF— exchange, 
virtual top quark contributions and contributions from other heavy particles such as W, Z 
bosons, charged Higgs particles, super symmetric particles in the supersymmetric extensions of 
the SM and other heavy objects in numerous extensions of this model. Consequently Ci{fi) 
depend generally on mt and also on the masses of new particles if extensions of the SM are 
considered. This dependence can be found by evaluating so-called box and penguin diagrams 
with full W, Z, top quark and new particles exchanges and properly including short distance 
QCD effects. The latter govern the /i-dependence of Ci{fj,). In models in which GIM mechanism 
[24] is absent, also tree diagrams can contribute to flavour changing neutral current (FCNC) 
processes. The point is that a given Cj receives generally contributions from all these three 
classes of diagrams. 

The value of /x can be chosen arbitrarily but the final result must be ;U-independent. Therefore 
the //-dependence of Cj(//) has to cancel the //-dependence of (Qj(/i)). In other words as far 
as heavy mass independent terms are concerned it is a matter of choice what exactly belongs 
to Cj(/i) and what to {Qi{fj,)). This cancellation of the //-dependence involves generally several 
terms in the expansion in (|2.2p . The coefficients Ci{^) depend also on the renormalization 
scheme. This scheme dependence must also be canceled by the one of (Qi(//)) so that the 
physical amplitudes are renormalization scheme independent. Again, as in the case of the //- 
dependence, the cancellation of the renormalization scheme dependence involves generally several 
terms in the expansion (j2.2p . One of the type of scheme dependences is the manner in which 
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75 is defined in D = 4 — 2e dimensions implying for instance the three schemes NDR, HV and 
DRED mentioned earher. 

Although fi is in principle arbitrary, it is customary to choose fi to be of the order of the 
mass of the decaying hadron. This is 0{m\y) and 0{mc) for B decays and D decays respectively. 
In the case of K decays the typical choice is /i = 0(1 — 2 GeV) instead of 0{mK), which is 
much too low for any perturbative calculation of the couplings Cj. Now due to the fact that 
/i <C M\Y^z, rnti large logarithms IuMw/m compensate in the evaluation of Ci{fi) the smallness 
of the QCD coupling constant Og and terms a"(lnMw//u)", {In M-\^ / fiy^'^ etc. have to be 
resummed to all orders in Og before a reliable result for Ci can be obtained. This can be done 
very efficiently by means of the renormalization group methods. The resulting renormalization 
group improved perturbative expansion for Ci{fi) in terms of the effective coupling constant 
OsifJ-) does not involve large logarithms and is more reliable. The related technical issues are 
discussed in detail in [T5] and [IB] and we will recall here only those that are essential for our 
presentation. 

All this looks rather formal but in fact should be familiar. Indeed, in the simplest case of 
the /3-decay, Heff takes the familiar form 

niff = ^ cosee[n7M(l - 75)d ® 67^(1 - 75)z^e] , (2.4) 

where Vud has been expressed in terms of the Cabibbo angle [19]. In this particular case the 
Wilson coefficient is equal unity and the local operator, the object between the square brackets, 
is given by a product of two V — A currents. Equation (|2.4p represents the Fermi theory for 
/3-decays as formulated by Sudarshan and Marshak [25] and Feynman and Gell-Mann |26j more 
than fifty years ago, except that in ()2.4p the quark language has been used and following Cabibbo 
a small departure of Vud from unity has been incorporated. In this context the basic formula 
()2.ip can be regarded as a generalization of the Fermi Theory to include all known quarks and 
leptons as well as their strong and electroweak interactions as summarized by the SM. 

Due to the interplay of electroweak and strong interactions the structure of the local operators 
is much richer than in the case of the /3-decay. They can be classified with respect to Lorentz 
structure, Dirac structure, the colour structure and the type of quarks and leptons relevant for 
a given decay. We will now list all the operators who's Wilson coefficients will be mentioned in 
subsequent sections. 

2.2 Local Operators in the SM 

We give below first a list of operators that play the role in weak B decays. Typical diagrams 
in the full theory from which these operators originate are shown in Fig. [TJ The cross in the 
diagram Id indicates that magnetic penguins originate from the mass-term on the external line 
in the usual QCD or QED penguin diagrams. The operators relevant for K decays are discussed 
subsequently. 
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u,c,t 




Figure 1: Typical Tree, Penguin and Box Diagrams in the SM. 
2.2.1 Nonleptonic Operators 

Of particular interest are the operators involving quarks only. In the case of the AS = 1 
transitions the relevant set of operators is given as follows: 
Current— Current (Fig. la): 



Ql = {Cab/3)v-A isi3Ca)v-A Q2 = {cb)v-A {sc)v-A 

QCD-Penguins (Fig. lb): 

Qs = {sb)v-A XI i^l)v-A Qa = {Sahp)v-A ^ {q^QaW-A 
q=u,d,s,c,b q=u,d,s,c,b 



(2.5) 



(2.6) 



Qb = {sb)v-A X {qq)v+A Qe = (sa6^)y_A X {qpqaW+A (2.7) 

q=u,d,s,c,b q=u,d,s,c,b 

Electroweak Penguins (Fig. Ic): 

<37 = 2 {sb)v-A XI ^1 (^qW+A Q»= 2 (^abfi)v-A X eqiQ/SQa)v+A (2.8) 
q=u,d,s,c,b q=u,d,s,c,b 



= 2 ^ eq{qq)v-A Qio = isabi3)v-A ^ eg {qpqa)v-A (2.9) 

q=u,d,s,c,b q=u,d,s,c,b 

Here, a, (3 denote colours and denotes the electric quark charges reflecting the electroweak 
origin of Q7, . . . , Qio- Finally, {cb)v-A = Cq, 7^(1 - 75)60,. 

These operators play a crucial role in non-leptonic decay of Bg and mesons and have 
through mixing under renormalization also an impact on other processes as we will see below. In 
this context let me also make one useful remark. In the literature operators in B physics appear 
sometimes with {hs)v-A or {sb)v-A, dependently whether B^ or B^ are studied, respectively. 
When using their Wilson coefficients given in the literature it is crucial to remember that they 
are complex conjugates of each other. This distinction is crucial for obtaining correct CP 
asymmetries. 

For non-leptonic K decays the flavours have to be changed appropriately. Explicit expressions 
can be found in [151 IS]. In particular the analogs of Qi and Q2 govern the A/ = 1/2 rule 
in Kl — )• TTvr decays, while the corresponding QCD penguins and electroweak penguins enter 
directly the ratio e'/e. 

2.2.2 Dipole Operators 

In the case of — Xs^ decay and B — )• Xgl^l^ decays as well as corresponding exclusive decays 
the crucial role is played by 

Magnetic Penguins (Fig. Id): 

= ^m^,So,cJ^'^{l + j5)baFf,, Qsg = ^m.So^a^" {I + i5)T^pbpG% (2.10) 

Again, when using the results in the literature care must be taken whether 6 or 6 is present in 
the operator and what are the factors multiplying the Dirac structures. The operator Qsg can 
also be relevant in nonleptonic decays. Also magnetic penguins with (1 + 75) replaced by (1 — 75) 
are present but they are suppressed within the SM with respect to the operators in ()2.10p by 

2.2.3 AF = 2 Operators 

In the case of — mixing and 5° - 5° mixing the relevant operators within the SM are 
A5 = 2 and AS = 2 Operators (Fig. le): 

Q(AS = 2) = {sd)v-A{sd)v-A Q{AB = 2) = (bd)v^A(bd)v-A (2.11) 
For Bg — Bg mixing one has to replace d by s in the last operator. 
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2.2.4 Semileptonic Operators 

In the case of i? — t- XJ^l^ also the following operators originating in Fig. If on top of magnetic 
penguins contribute 

Qgv = isb)v-AiP'fJ')v QlOA = isb)v-A{p'fJ')A- (2.12) 

Changing appropriately flavours one obtains the corresponding operators relevant for B — t- 
Xdl+r and Kl ir^l+r 

The rare decays B — t- Xgi'i^, B — t- K*vi>, B — t- Kvu and Bg — s- Ji^ are governed by 

Qvv{B) = {sb)v-A{vv)v-A QtJ.p.iB) = {sb)v-Aifif^)v-A ■ (2-13) 

The rare decays K — )• vrz^P and Kl — )• flfi are governed on the other hand by 

Qvv{K) = {sd)v~A{vy)v-A Qtifi{K) = {sd)v^A{Ulj)v~A ■ (2-14) 

2.3 Local Operators in Extensions of the SM 

New physics (NP) can generate new operators. Typically new operators are generated through 
the presence of right-handed (RH) currents and scalar currents with the latter strongly sup- 
pressed within the SM. New gauge bosons and scalar exchanges are at the origin of these op- 
erators that can have important impact on phenomenology. Below we give examples of new 
operators being aware that this list is incomplete. Much more extensive discussion can be found 
in [271 [28]. 



2.3.1 AF = 2 Non-leptonic Operators 

For definiteness, we shall consider here operators responsible for the K^-K^ mixing. There are 
8 such operators of dimension 6. They can be split into 5 separate sectors, according to the 
chirality of the quark fields they contain. The operators belonging to the first three sectors 
(VLL, LR and SLL) read [29] (our competition in Rome [30] uses a different basis): 

QY^ = {-s''l^.PLd''){s^rPRd^). 
qLR = {-s-PLd^){-s^PRd^), 

gSLL ^ (s°Pid")(s^PLdO, 

qSLL ^ {ra^,PLd^){-s^a>'''PLd^), (2.15) 

where a^y = ^[7^,7^/] and Pl,r = |(1 =F 75)- The operators belonging to the two remaining 
sectors (VRR and SRR) are obtained from Q^^^ and Qf^^ by interchanging P^ and Pr. In the 
SM only the operator QJ'^^ = Q{AS = 2)/4 is present. 
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2.3.2 AF = 1 Non-leptonic Current-Current Operators 

In the present section, we list the current-current four-quark AF = 1 operators. For this 
purpose, we choose the operators in such a manner that all the four flavours they contain are 
different: s, d, u, c. In such a case, the only possible diagrams are the current-current ones. 
Penguin diagrams are discussed subsequently. 

Twenty linearly independent operators can be built out of four different quark fields. They 
can be split into 8 separate sectors, between which there is no mixing. The operators belonging 
to the first four sectors (VLL, VLR, SLR and SLL) read [29] 



Q^" = 






g^^ = 

^VLR 


{s''l^PLd'^){vP^^PRcf') = Qv.Vn, 




gf ^ = 

^2 — 


{s^PLdP){uPpRn = Qlr, 
is''PLd''){u^PRC^) = Qlr, 




gf ^ = 
gi^^ = 

V4 — 


{s^PLdP){uPPLn = QlL, 

{s''PLd''){u^PLC^) = Qll, 
{s''a^,PLd^){u''a^''PLn = Qt^t,, 
is''a^,PLd''){u^a^'''PLC^) = Qt,t„ 


(2.16) 



where on the r.h.s. we have shown the notation of the Rome group |30j . 

The operators belonging to the four remaining sectors (VRR, VRL, SRL and SRR) are 
obtained from the above by interchanging Pl and Pr. Obviously, it is sufficient to calculate the 
anomalous dimensions (ADMs) only for the VLL, VLR, SLR and SLL sectors. The "mirror" 
operators in the VRR, VRL, SRL and SRR sectors will have exactly the same properties under 
QCD renormalization. On the other hand their Wilson coefficients being governed by weak 
interactions can be different. In the SM only the operators QY^^ and Q^^^ are present. 

2.3.3 A.F = 1 Non-leptonic Penguin Operators 

The operators in (j2.15p and (j2.16p do not constitute the full set of six-dimensional four quark 
operators contributing to AF = 1 processes. In addition to QCD penguins and electroweak 
penguins of the SM there are other penguin operators. In our paper |29] we have therefore 
generalized our analysis of two-loop anomalous dimensions to the full set of AF = 1 four-quark 
operators. These results are much less known but should be useful in the extensions of the SM 
one day. The list of these operators can be found in [25] . 
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2.3.4 Dipole Operators 

In the presence of right-handed (RH) currents, mediated for instance by a very heavy Wr in 
left-right symmetric models the magnetic penguins 

^77 = ^rubSaa^^-'il - j5)baFf,, QsG = ^m6S,(7'^-(l - 75)T^0bpG% (2.17) 

could be important. 

2.3.5 A.F = 1 Semi-leptonic Operators 

Concerning the semi-leptonic operators in the extensions of the SM the typical examples of 
operators related to the presence of RH currents are 



Qqv = {sb)v+A{tJ'fJ')v QioA = {sb)v+A{tJ'fJ')A- (2.18) 

Quv{B) = {sb)v+A{i>i^)v~A Q^l^i{B) = (s6)v/+A(w)y-A • (2.19) 

Quv{K) = {sd)v+A{vv)v-A QfifiiK) = {sd)v+A{fJ'tJ')v-A ■ (2.20) 

If scalar currents resulting from scalar exchanges like the heavy Higgs in the 2HDM models 
or sparticles in the MSSM are present, scalar operators enter the game. The most prominent 
are the ones that govern the Bg — )• decay in 2HDMs and the MSSM at large tan/3: 

Qs = {sPLh){m) Qp = {sPLb){jll^y.). (2.21) 

Qs = {sPRb){m) Qp = isPRb)ifl-f5f^). (2.22) 

2.4 Wilson CoefRcients 
2.4.1 General Structure 

The main objects of interest in this review are the QCD and electroweak corrections to the 
Wilson coefficients of the operators listed above. Once these coefficients have been calculated 
at a high energy scale like Mw §], the renormalization group methods allow to calculate them 
at low energy scales at which the matrix elements are evaluated by means of non-perturbative 
methods. Denoting this lower scale simply by fi the general expression for Ci{fi) is given by: 

C{fi) = U{fi,Mw)C{Mw) , (2.23) 



^Ci{Mw) are often called matching conditions as they are found through matching of the full theory with 
heavy fields as dynamical degrees of freedom to the effective theory where only light fields are dynamical. 
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where C is a column vector built out of Cj's. C{Mw) are the initial conditions which depend 
on the short distance physics at high energy scales. In particular they depend on mt and the 
masses and couplings of new heavy particles in the extensions of the SM. We set the high energy 
scale at Mw, but other choices are clearly possible. U{fi, Mw), the evolution matrix from M]y 
down to n, is given as follows: 



U{iJ,,Mw) = Tgexp 



'''' dg'f^ 



(2.24) 



with g denoting the QCD effective coupling constant and Tg an ordering operation defined in 
[16j . 13 {g) governs the evolution of g and 7 is the anomalous dimension matrix of the oper- 
ators involved. The structure of this equation makes it clear that the renormalization group 
approach goes beyond the usual perturbation theory. Indeed U{^, Mw) sums automatically 
large logarithms log M^/fi which appear for /i <C M^r. In the so-called leading logarithmic 
approximation (LO) terms (g^ log Mi4///i)" are summed. The next-to-leading logarithmic cor- 
rection (NLO) to this result involves summation of terms (g^)"(logMi4///x)"~^ and so on. This 
hierarchic structure gives the renormalization group improved perturbation theory. 

As an example let us consider only QCD effects and the case of a single operator so that 
({2:23]) reduces to 



Ciii) = U{fi,Mw)C{Mw) 

with C(//) denoting the coefficient of the operator in question. 

Keeping the first three terms in the expansions of 7(5) and I3{g) in powers of g: 



(2.25) 



47r 



\47r/ 



in 



(2.26) 



/3(5) = -/3o 



9 



/3i 



9 



9 



'I67r2 ^'(167r2)2 ^^(167r2)3 
and inserting these expansions into (j2.24p gives: 



(2.27) 



U{fi,Mw) 



H : Ji + 



in 



in 



J2 



w 



as{Mw) , 
^ i^-^'- 



as{Mw)\ , 2^ 

[J2-JI) 



in 



(2.28) 



where 



P 



J2 



2/3o' ' 2/3o 



P 

2A 



^(2) 

Wo 



(2.29) 
(2.30) 
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General formulae for the evolution matrix U{fi, M\\r) in the case of operator mixing and valid also 
for electroweak effects at the NLO level can be found in [15]. The corresponding NNLO formulae 
are rather complicated and given in [31] • The leading logarithmic approximation corresponds 
to setting = J2 = in ^(TM) . In the NLO only J2 = and the last term in ^^TM) has to be 
removed. 

The coefficients /3j are given as follows 

33-2/ „ 306 - 38/ 
/3o = ^ /3i = (2.31) 



2 18 " 54 
where / is the number of quark flavours. 
The expansion for C(Mw) is given by 

C(Mw) = Co + ^^C, + (^^) ' (2.33) 

where Co, Ci and C2 depend generally on mt, M^i/, the masses of the new particles and the new 
parameters in the extentions of the SM. It should be stressed that the renormalization scheme 
dependence of Ci and C2 is canceled by the one of Ji and J2 in the last square bracket in (I2.28P 
although at the NNLO level this cancellation is rather involved. The scheme dependence of Ji 
and J2 in the first square bracket in (j2.28p is canceled by the scheme dependence of {Q{fi)). 
The power P is scheme independent. The methods for the calculation of U{fi, Mw) and the 
discussion of the cancellation of the and renormalization scheme dependences are presented 
in detail in [16] and in the original papers where such calculations have been done. 

When talking about the ^-dependence one should distinguish two types of dependences. 
First we have the // dependence related to the renormalization of operators and present in the 
evolution matrix. This dependence arises in the presence of non- vanishing anomalous dimensions 
of the operators responsible for weak decays. But in the coefficients Ci{Mw) in (j2.33p also heavy 
quark masses are present that are running masses with their scale dependence governed by the 
anomalous dimension of the mass operator 

7„M = 7!»'£ + 7«(^)%7S'(£)' P'^^) 
with the coefficients 'ym given as follows [32l [331 IS] 



^(2) 

im 



/2216 160 „ 140 „2 



(2.36) 
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where C(3) ~ 1.202057. These results are vahd in the MS scheme. The four-loop contribution 
has been found in [35] but it is too complicated to be presented here. The most recent values of 
running quark and lepton masses can be found in [36] . 

If only the leading term Co is present the choice of the /i for the masses matters. This 
unphysical scale dependence is cancelled by the non- leading terms in (|2.33|) . A detailed discussion 
of this issue can be found in [16] and we will return to it briefly below. 

For later purposes it will be useful to generalize the formula (j2.26p to include mixing between 
operators and 0{a) effects, where a is the QED coupling constant. This formula is relevant 
whenever also electroweak effects are present and electroweak penguin operators contribute. At 
the NNLO level in QCD but to leading order in a one has: 

with 7^'^^ , 7s^^ and 7^^'* being anomalous dimension matrices that are generalizations of the cor- 
responding coefficients in (j2.26p to include mixing among operators under QCD renormalization. 
liO{a) effects are included in the coefficients at scales O(-Mw), the anomalous dimension matrix 
must also include 0{a) contributions which are represented by 76*^^ ji]) and 7se^ at LO, NLO 
and NNLO, respectively. 

The corresponding generalization of the Wilson coefficients in (j2.33p takes the form 

C(Mw) = Co + ^^C, + [^-^^) C2 + ^Ci + ^^^C, (2.38) 

where now the coefficients are column vectors and the evolution in (I2.25P generalizes for a = 
to the one in ()2.23p . For a ^ similar formulae exist. We will give later the directions to papers, 
where explicit expressions for all these objects can be found. 

2.4.2 Renormalization Scheme Dependence 

As already stated above, beyond LO various quantities like Wilson coefficients and the anomalous 
dimensions depend on the renormalization scheme for operators. This dependence arises because 
the renormalization prescription involves an arbitrariness in the finite parts to be subtracted 
along with the ultraviolet singularities. Two different schemes are then related by a finite 
renormalization. 

I have discussed this issue in detail in [16J, in particular in Section 6.7 of these lectures. Here 
I just want to recall one NLO formula to which I will refer from time to time. It is a relation 
between anomalous dimension matrices in two different renormalization schemes: 

^(0)/ ^ ^(0) ^(1)/ ^ ^(1) ^ [^^^ ^(0)] ^ 2/3o Ar, (2.39) 

where the prime distinguishes the two schemes and Af is a shift at 0{as) in the matrix elements 
of operators calculated in these two renormalization schemes: 

(Q)' = (1 + ^Af)(g), C' = (1 - ^(Af)^)C. (2.40) 
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2.5 Inclusive Decays 

So far I have discussed only exclusive decays. It turns out that in the case of inclusive decays 
of heavy mesons, hke S-mesons, things turn out to be easier. In an inclusive decay one sums 
over all (or over a special class) of accessible final states and eventually one can show that the 
resulting branching ratio can be calculated in the expansion in inverse powers of mb with the 
leading term described by the spectator model in which the S-meson decay is modelled by the 
decay of the 6-quark. Very schematically one has then for the decay rate 



This formula is known under the name of the Heavy Quark Expansion (HQE) [S^ . Pedagogical 
reviews on this topic and heavy quark effective theories can be found in |371 [381 ES HQl H]. 

Since the leading term in this expansion represents the decay of the quark, it can be calculated 
in perturbation theory or more correctly in the renormalization group improved perturbation 
theory. It should be realized that also here the basic starting point is the effective Hamiltonian 
(j2.1|) and that the knowledge of the couplings Cj (/i) is essential for the evaluation of the leading 
term in (j2.4ip . But there is an important difference relative to the exclusive case: the matrix 
elements of the operators Qi can be "effectively" evaluated in perturbation theory. This means, 
in particular, that their and renormalization scheme dependences can be evaluated and the 
cancellation of these dependences by those present in Ci{fi) can be investigated. 

Clearly in order to complete the evaluation of T{B — t- X) also the remaining terms in (I2.4ip 
have to be considered. These terms are of a non-perturbative origin, but fortunately they 
are suppressed often by two powers of rrib. They have been studied by several authors in the 
literature with the result that they affect various branching ratios by less than 10% and often 
by only a few percent. Consequently the inclusive decays give generally more precise theoretical 
predictions at present than the exclusive decays. On the other hand their measurements are 
harder. There are of course some important theoretical issues related to the validity of HQE 
in (|2.41|) which appear in the literature under the name of quark-hadron duality but I will not 
discuss them here. 

The very rough appearance of the second term on the r.h.s of (j2.4ip totally underrepresents 
the efforts which have been made over many years to calculate these contributions. But as I was 
not involved in these efforts and they contain some non-perturbative aspects, I will not discuss 
them here. On the other hand we will summarize the status of the first term in ()2.4ip in some 
detail in Sections 6 and 8. 




1 



(2.41) 



■h 
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2.6 The Structure and the Status of the NLO and NNLO Corrections 



2.6.1 General Comments 

As we will see in the following sections during the last two decades the NLO corrections to Ci {fi) 
have been calculated within the SM for the most important and interesting decays. Also several 
NNLO calculations have been performed. In tables [2][6] we give references to all NLO and NNLO 
calculations within the SM done until the end of February 2011 that deal with the processes 
discussed by us. While these calculations improved considerably the precision of theoretical 
predictions in weak decays and can be considered as an important progress in this field, the 
pioneering LO calculations for current-current operators [l2l |33], penguin operators [lH H5] . 
A5 = 2 operators [l6] and rare K decays [47J should not be forgotten. 

2.6.2 Different Classes of QCD Corrections 

The structure of QCD corrections to various decays depends on the decay considered. In par- 
ticular the expansion in as can vary from decay to decay. Moreover even within a given decay 
the structure of QCD corrections to internal charm and top contributions differ from each other. 
Let us then classify various cases beginning with the simplest ones and systematically increasing 
the complexity. 
Class 1 

The simplest situation arises when there is only one contributing operator with a vanishing 
anomalous dimension and in addition there is no mixing of this operator with operators which 
have non-vanishing anomalous dimensions. Moreover the diagrams in the full theory from which 
this operator was born have only heavy internal particles like W^, the top quark and generally 
heavy particle exchanges in the extensions of the SM. 

This is the case of the operators Quu{B), Q^p_{B) in (j2.13p and Quu{B), Qfj,ii{B) in (|2.20p 
contributing to rare decays B — t- XsUi), B — t- K*vi>^ B — t- KvT> and Bg — s- /x^/i^. If the 
operators QuuiK) and Q^fi{K) in (j2.14p and Quu{K) and Q^^{K) in (|2.20p originate in the 
internal top quark contributions and other heavy particle contributions to rare decays K — )• ttuv' 
and Kl — t- h'^ then also these contributions belong to this class. The case of internal charm 
quark contributions is classified separately below. 

Denoting a given loop function in the absence of QCD corrections by Fi[x), the decay 
amplitudes in this case have the following perturbative expansion in 

Ai = Fi{x) + 0{as) + 0{al) (2.42) 

with as evaluated at the high scale where the operator has been generated. We drop the Lorentz 
structure for simplicity. Therefore 0{as) corrections are generally small and no large logarithm 
related to operator renormalization is present in them due to the vanishing of the anomalous 
dimension of the contributing operator. However x = rn^/My^ depends on a scale jit present 
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in mt{^t) with similar comments applying to other coloured heavy particles present beyond the 
SM. The corresponding logarithm involving this scale and present in the 0{as) correction in 
(j2.42p cancels this scale dependence present in the leading term Fi (x) so that up to higher order 
corrections Ai is independent of fif On the other hand the size of the 0{as) correction in (j2.42p 
clearly depends on the chosen fit- It turns out that it is useful to set fit = rnt{mt). Then the 
result can be summarized by 

^1 =Fi(x)r/QCD (2.43) 

with ryqcD close to unity and practically independent of the measured top quark mass. For other 
choices of fit the factor ?7qcd can defer significantly from unity but then also the numerical value 
of Fi[x) is different so that Ai remains the same up to higher order corrections. This removal 
of order 10% dependence on fit in the LO formulae for rare K and B decays was the basic 
motivation for the calculations in \4S\ H9] . The QCD calculations in this class are described in 
Section 5 and the relevant references are collected in Table HI 
Class 2 

This class is constructed from Class 1 by giving the operator an anomalous dimension but 
still requiring that it does not mix with other operators and all particles in loops generating this 
operator are heavy. This is the case of AS" = 2 and /S.B = 2 operators in the SM when only top 
quark contributions in the box diagrams are considered. In this case we have 

A2 = F2{x) * Iqcd + Oias) + O(a^), (2.44) 

where the funny Iqcd represents the leading RG factor like the one involving P in (|2.28|) . Now 
the 0{as) and higher order terms involve not only Mw but also low energy scale fi at the end of 
the RG evolution as seen explicitly in (j2.28p . Moreover the 0{as{Mw)) correction involves now 
two logarithms multiplied by two different anomalous dimensions, one anomalous dimension of 
the mass operator related to the fit dependence present already in Class 1 and the second present 
in P involving the anomalous dimension of Q{AF = 2) in (j2.1ip . The latter logarithm cancels 
the fiw dependence present in the funny factor Iqcd in (|2.44p so that A2 does not depend on the 
precise value of the scale at which the Wilson coefficients are defined. Again one can summarize 
the result schematically by 

^2 = i^2(a;)?7QCD. (2.45) 

However, this time ??qcd can depart significantly from unity as summation of large logarithms in 
the process of RG evolution is involved. ?7qcd depends as seen in (j2.28p on the lower scale fi and 
this dependence cancels the one present in the hadronic matrix elements. This later dependence 
in 77QCD is often factored out so that the known factors 772 = r]tt in ex and rfs in ^ — ^ 
mixing are //-independent and this also applies to the Bi factors that up to factors involving 
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weak decays constants represent hadronic matrix elements. Explicit expressions are given in 
Section 4. 

This discussion applies also to the AF = 2 operators in (j2.15p except that now mixing under 
renormalization between operators and and similarly between Qf^^ and Q^^^ takes 
place. Explicit formulae for this case can be found in [50]. 

In some extensions of the SM FCNC operators are generated already at tree level but also 
in this case analogous discussion can be made. 

The QCD calculations in this class are described in Section 4 and the relevant references are 
collected in Table El 

Class 3 

We next consider QCD corrections to charm contributions to — )■ ir^i^v' and i^L — ^ ^^/^^ 
decays. The corresponding operators are those of Class 1 but now a light charm quark mass is 
present in the loop. Consequently even without QCD corrections a large logarithm Innic/Mw 
is present and on the way to low scales bilocal operators enter the game. They undergo a rather 
complicated renormalization [5Tj. Only when the charm is integrated out we obtain the local 
operators Qi,p{K) and Q^ji{K) in (j2.14p . From the point of view of the renormalization group 
analysis, the expansion in as takes in this case the following form 

A3 = 0(— ) + + (2.46) 

as 

Thus the NLO corrections to the charm part of the amplitudes for — t- tt^uv and Ki^ — )• 
amount to the 0{1) term (51], while the NNLO corrections to the 0{as). Still to get this term 
three-loop diagrams have to be evaluated [311 [52]. The LO term has been calculated in [37]. 

The QCD calculations in this class are described in Section 5 and the relevant references are 
collected in Table HI 

Class 4 

We next consider the set of AF = 1 current-current operators in (j2.5p that both have 
nonvanishing anomalous dimensions and mix under renormalization. After diagonalization of 
this system one gets two operators Q± in (13. 8p which evolve without mixing. The situation for 
each operator is then similar to Class 2 except that no box diagrams involving heavy particles 
have to be evaluated. Consequently the expansion is as follows 

A4 = Iqcd + O (as) + 0{al), (2.47) 

where as terms are evaluated both at Myy and the low scale ^ according to the evolution in 
()2.28p . The leading term is again the one involving P in ()2.28p . 

The QCD calculations in this class are described in Section 3 and the relevant references are 
collected in Table [2j 
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Class 5 

We next consider QCD penguin and electroweak penguin operators in (I2.6p - (l2.9p contributing 
to non-leptonic decays. These operators mix under QCD and QED renormalization and evidently 
the QCD penguin and electroweak penguin diagrams in the full theory are 0{as) and 0{a), 
respectively. The Wilson coefficients of the corresponding operators after the top quark and 
have been integrated out are also of the same order respectively. This mismatch of powers in ag 
can be overcome in the process of renormalization group analysis by multiplying the electroweak 
operators by l/a^ and compensating this rescaling by multiplying their Wilson coefficients by 
as- As a is from the point of view of QCD a fixed number, the RG can now be performed as in 
class 2 except for the following changes. Qi and Q2 operators have to be included as they mix 
into Q3 — Qio operators affecting their QCD evoloution. Thus we deal with 10 x 10 anomalous 
dimension matrices but because of the presence of electroweak penguins also terms 0{aas) at 
NLO order have to be considered and 0{aa'^) at NNLO in addition to the usual 0{a'^) and (af ) 
terms, respectively. See (j2.37p . 

We observe that now the formulae are a bit more involved but what is more important are 
the following facts which apply for instance to the evaluation of the ratio e'/e in Kl — )• vrvr: 

• At LO there is no top quark mass dependence nor any heavy particle mass dependence 
from penguin diagrams. 

• At NLO these mass dependences enter for the first time. 

• The NLO matching conditions for electroweak penguin operators do not involve QCD 
corrections to box and penguin diagrams and consequently the renormalization scale de- 
pendence in the top quark mass in these processes is not negligible at the NLO. In order 
to reduce this unphysical dependence, QCD corrections to the relevant box and penguin 
diagrams have to be computed. In the renormalization group improved perturbation the- 
ory these corrections are a part of NNLO corrections. In [S3] such corrections have been 
computed for e'/e. Their inclusion allowed to reduce renormalization scheme dependence 
present in the electroweak penguin sector. 

The QCD calculations in this class are described in Section 3 and the relevant references are 
collected in Table [2j 
Class 6 

We next come to magnetic penguin operators in (|2.10p restricting the discussion to the 
B —7- Xsj decay. These two operators mix under renormalization with each other and are also 
influenced by the mixing with the current-current operators and QCD penguin operators. Thus 
we deal here with a 8 x 8 anomalous dimension matrix. 

I will report on the heroic efforts to calculate the QCD corrections to the B — )■ Xg^ rate in 
Section 6. Let me here only write down symbolically the general structure of the amplitudes in 
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this class: 



^6 = Fq{x) * Iqcd + Iqcd + 0{as) + ©(a^), 



(2.48) 



with the first two terms representing LO and the 0{as) and 0{a'^) terms NLO and NNLO 
corrections, respectively. The Iqcd term represents the mixing of magnetic penguin operators 
with the current-current and QCD-penguin operators. In fact this term is responsible for a 
strong enhancement of the B — )• decay rate. 

The QCD calculations of — t- Xg^ decay rate are described in Section 6 and the relevant 
references are collected in Table [5j 



In rare decays Kl — >• tt^I^I^ and B — t- XJ^l^ semileptonic operators in ()2.12p are involved. 
In addition to these operators also current-current non-leptonic operators (j2.5p and QCD pen- 
guin operators in (j2.6p and (|2.7p have to be taken into account. The electroweak penguins and 
magnetic penguins turn out to be irrelevant in — )• 'k^I~^1~ but the magnetic penguins have 
to be included in i? — )• Xsl'^l~ ■ 

The new feature is that although the operator Q^y has no anomalous dimension by itself it 
mixes with the operators Qi — Qq and consequently a 7 x 7 anomalous dimension matrix has to 
be considered. The resulting structure of the coefficient Cgy is then 

C9v = 0(— ) + F(x) + iQCD + 0(a.). (2.49) 

Consequently the structure is similar to Class 3 but at the NLO level in this case heavy quark 
mass dependence enters. 

The operator Qiqa has no anomalous dimension and similarly to operators in (|2.13p and 
(j2.14p does not mix with anybody. Therefore its Wilson coefficient has the structure in (j2.42p . 
Evidently this operator enters the amplitude for — )• 7r^l~^l~ at the NLO level. 

In the case of -B — )• Xsl~^l~ the situation is complicated by the presence of magnetic operators 
but the structure of the corresponding Wilson coefficients Cgy and Ciqa is the same. The 0{as) 
corrections to the penguin and box diagram relevant for the NNLO evaluation of -B — t- Xgl^l^ 
rate have been calculated in [55]. 

The QCD calculations of Kl 7r°Z+/- and B Xsl+l~ are described in Sections 7 and 8 
respectively and the relevant references are collected in Table [6l 



Finally I give the structure of QCD corrections for the charm-charm (rji = rjcc) and charm- 
top (r/3 = rjct) contributions to AS = 2 Hamiltonian. The structure of these corrections differs 
from r]2 = rju discussed in Class 2 and also from each other but I think it is instructive to put 
them together in order to see the difference. We have 



Class 7 



Class 8 



Vi = {asf+ (Iqcd + 0(a.) + 0(a' 



)) 



(2.50) 
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m = (as) 




Iqcd + Iqcd + 0{as) 



) 



(2.51) 



The references to QCD calculations of iji and 1^3 are collected in Table [3] where also several 
remarks on these calculations will be made. 

2.6.3 Two-Loop Anomalous Dimensions Beyond the SM 

In the extentions of the SM new operators are present. The two loop anomalous dimensions 
for the AF = 2 and AF = 1 four-quark dimension-six operators listed in (I2.15P and (12.161) 
have been computed in \30\ I29j . In [29] also the remaining anomalous dimensions of AF = 1 
four-quark operators can be found. 

2.6.4 Tvi^o-Loop Electroweak Corrections 

In order to reduce scheme and scale dependences related to the definition of electroweak pa- 
rameters like sin^ 9w, and oqed, two-loop electroweak contributions to rare decays have to be 
computed. For — )• tt^uu, B^.s and B — >■ Xsiyv they can be found in [56^ I57j . for 

Bd.s - Bd,s mixing in [58j and for B Xs^i in [59l EQI EH ESI [63] • 

2.6.5 NLO QCD Calculations Beyond the SM 

There exist also a number of partial or complete NLO QCD calculations within the Two-Higgs- 
Doublet Model (2HDM) and the MSSM. In the case of the Two-Higgs-Doublet Model such 
calculations for B^^ — B^^ mixing, B — )• X^^ and B — )• Xsl~^l~ can be found in [63|, [55 [ [65 | [66 l 
[671 [68] and [69], respectively In 2HDM also NNLO QCD Corrections to B ^ X^'j have been 
recently calculated [70] . 

The corresponding NLO calculations for B'^ ^ — B^ ^ and B — )■ Xg'y in the MSSM can be found 
in[nill2l[73l|71]and[5a[7a|76llI71[78l|79], respectively The paper [55] gives also the results 
for B —7- Xsgluon. In fact Bobeth, Misiak and Urban [55] present rather general formulae for 
Wilson coefficients relevant for B — t- Xg^ and B ^ Xg gluon evaluated at high scale (matching 
conditions) at the LO and NLO level that can be used for other extensions of the SM. 

Finally, I would like to mention two calculations in which I took part: NLO QCD corrections 
to rare K and B decays in the MSSM at low tan (3 [8OJ and NNLO QCD corrections to S — )• 
Xsl+l- in the MSSM [ST]. 

2.6.6 Penguin-Box Expansion 

The rare and CP violating decays of K and B mesons as well as ek, e'/e and - -8° mixmgs 
are governed in the SM by various penguin and box diagrams with internal top quark and 
charm quark exchanges. Some examples are shown in Fig. [TJ Evaluating these diagrams one 
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finds a set of basic universal (process independent) mt-dependent functions Fr{xt) [82l[T6] where 
/M^. Explicit expressions for these functions can be found in [16j. 
It is useful to express the OPE formula ()2.2p directly in terms of the functions Fr{xt) |83j : 



AiM ^F) = Po{M ^ F) + Y^ Pr{M F) Fr{xt), (2.52) 

r 

where the sum runs over all possible functions contributing to a given amplitude. Pq summarizes 
contributions stemming from internal charm quark. In the OPE formula (|2.2|) . the functions 
Fr{xt) are hidden in the initial conditions for Ci{fi) represented by (7(M-w) in (|2.23p . 

The coefficients Pq and Pr are process dependent and include QCD corrections contained 
in the evolution matrix U{iJ,, My^). They depend also on hadronic matrix elements of local 
operators and the relevant CKM factors. An efficient and straightforward method for finding 
the coefficients Pr is presented in [83]. As the expansion in (j2.52p involves basic one-loop 
functions from penguin and box diagrams it was naturally given the name of the Penguin-Box 
Expansion (PBE). 

Generally, several basic functions contribute to a given decay, although decays exist which 
depend only on a single function. We have the following correspondence between the most 
interesting FCNC processes and the basic functions within models with constrained Minimal 
Flavour Violation: 

-mixing (ek) S{v) 
Bis - ^L-mixing (AM,,^) S{v) 
K — )- TTuv, B — ;> Xd^sT-'v X{v) 

Ki, fifi, Bd,s II y{v) 

vrOe+e- Y{v), Z{v), E{v) 

e', Nonleptonic = 1, AS* = 1 X{v), Y{v), Z{v), E{v) 

B ^ D'{v), E'{v) 

B ^Xs gluon E\v) 

B Xsl+l- Y{v), Z{v), E{v), D'{v), E'{v), 

where v denotes collectively the arguments of a given function with v = xt in the SM. In these 
models the operator structure of the SM remains intact and NP modifies only the basic functions. 
It should be mentioned that this correspondence is strictly valid in LO. At NLO in processes 
in which mixing between different operators is present new loop functions can contribute to a 
given process but these contributions are generally small. 

Originally PBE was designed to expose the mt-dependence of FCNC processes [83] which 
was hidden in the Wilson coefficients. In particular in the case of e' /e, where many of these 
functions enter, this turned out to be very useful. After the top quark mass has been measured 
precisely this role of PBE is less important. On the other hand, PBE is very well suited for 
the study of the extensions of the SM in which new particles are exchanged in the loops. If 
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there are no new local operators beyond those present in the SM the mere change is to modify 
the functions Fr{xt) which now acquire the dependence on the masses of new particles such as 
charged Higgs particles and supersymmetric particles. The process dependent coefficients Pq 
and Pr remain unchanged. The effects of new physics can then be seen transparently. Many 
examples of the applications of PBE can be found in the literature. In particular in the last 
decade we have used this method for the study of FCNC processes in the MSSM, a model with a 
universal extra dimension, littlest Higgs model (LH), littlest Higgs model with T parity (LHT), 
MSSM at low tan f3 and also the SM with four generations. In these papers compilations of the 
functions Fr in a given model can be found. A complete list of references can be found in [84J. 

One virtue of this method is a transparent study of the departure from MFV. In this frame- 
work, as discussed in detail in |85j . the basic loop functions are universal with respect to the 
system considered. Indeed as seen above the same function S{v) enters ek and Bd^s — B^^s 
mixings. Similarly the same function X{v) enters — )• tt'^vv' and B — )■ X^ gi^i' decays. Indeed 
in MFV the flavour dependence resides fully in the CKM matrix elements. Moreover, all these 
functions are real as the sole complex CP-violating phase resides in the CKM matrix and in 
flavour blind CP phases. In the presence of new sources of flavour and CP violation things are 
different: 

• The universality in question is broken and many MFV relations between various branching 
ratios are generally violated. 

• The basic functions become complex quantities leading often to new CP-violating effects 
not encountered in the CKM framework. 

All these new effects can be transparently seen in this framework. 

If new effective operators with different Dirac and colour structures are present, new functions 
multiplied by new coefficients Pr{M — )• F) contribute to (j2.52p . Still one is free to add these 
contributions to the functions of the SM, but then the dependence on the new Pr{M — )• F) is 
included in the basic functions. This is not always convenient if one wants to see explicitly the 
effects of new right-handed operators or scalar operators Qf ^, Qf^^ and Qf^^ given previously. 
Therefore in this case it is better to proceed as follows. 

We can start with ()2.2p but instead of evaluating it at the low energy scale we choose 
for fi the high energy scale to be called fj,H at which heavy particles are integrated out. Then 
absorbing Gp/^/2 and V^j^jy,^ in the Wilson coefficients Ci^fin) the amplitude for M — M mixing 
(M = K,Bd,Bs) is simply given by 



Here the sum runs over all the operators in (j2.15p . that is i = 1, 2 and a = VLL, VRR, LR, .... 
The matrix elements for — B^ mixing are for instance given as follows |50] 




(2.53) 



imiQtMlB', 



) = ImI,fIp-{b,), 



(2.54) 
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where the coefficients P^{Bd) collect compactly all RG effects from scales below fj,H as well as 
hadronic matrix elements obtained by lattice methods at low energy scales. Analytic formulae 
for these coefficients are given in [50j while the recent applications of this method can be found 
in [861 EZl EH]- As the Wilson coefficients Ci{fiH) depend directly on the loop functions and 
fundamental parameters of a given theory, this dependence can be exhibited if necessary. 

Again as in the case of PBE the virtue of using high energy scale rather than the low energy 
scale is that the coefficients P"(M) can be evaluated once for all if the hadronic matrix elements 
are known. 

The following points should be emphasized: 

• The expressions ()2.53p and ()2.54p are valid for any model with the model dependence 
entering only the Wilson coefficients Cf'(nH), which generally also depend on the meson 
system considered. In particular, they are valid both within and beyond the MFV frame- 
work. In MFV models CKM factors and Yukawa couplings define the flavour dependence 
of these coefficients, while in non-MFV models additional flavour structures are present in 

• The coefficients are model independent and include the renormalization group evolution 
from high scale fin down to low energy 0{fiK, tJ-s)- As the physics cannot depend on the 
renormalization scale ^h, the depend also on fin so that the scale dependence present 
in Pf is canceled by the one in Cf . Explicit formulae for fj^n dependence of can be 
found in |5Uj . It should be stressed that here we are talking about logarithmic dependence 
on fifj- The power-like dependence (such as . . . ) is present only in the C". 

• The P° depend however on the system considered as the hadronic matrix elements of 
the operators in (|2.53p relevant for instance for — mixing differ from the matrix 
elements of analogous operators relevant for B'^ ^ — ^ systems. Moreover whereas the 
RG evolution in the latter systems stops at fiB = 0{Mb), in the case of — system 
it is continued down to ~ 2 GeV, where the hadronic matrix elements are evaluated 
by lattice methods. 

After this rather heavy material we are ready to return to the main story of this paper. 

3 A/S = 1 and AS = 1 Non-Leptonic and Semi-Leptonic Decays 
3.1 Effective Hamiltonians 

The effective Hamiltonian for non-leptonic AS = 1 transitions is given in the SM as follows: 
G ^° 

T^eff (A5 = 1) = -^KTsKd E (^^(/^) + ^ ^«(^) (3.1) 
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with r = — V^*V^^/{V*gV^^). The operators Qi are the analogues of the ones given in (j2.5p - (|2.9p . 
In the case of AB = 1 transitions the flavours have to be changed appropriately and the effective 
hamiltonian is usually written as follows 



H,ff{AB = l) = ^ 



10,8G 



V2 



_ (3.2) 

i=3 

where 

A, = V*,Vgb (3.3) 

and 

Ql = {qabi3)v-Aisi3qa)v-A , Ql = {qaba)v-A{si3qi3)v-A ■ (3.4) 

In particular Ql = Qi in (j2.5p . 



3.2 Current-Current Operators 

Let me begin our NLO story with the first climb within the MNLC which Peter Weisz and 
myself started in December 1988: the calculation of two-loop anomalous dimensions of Qi and 
Q2 operators. It involves 28 two-loop diagrams shown in Fig. [3] and the corresponding counter 
terms. We decided to perform it in three schemes for 75: anticommuting (NDR), 't Hoof-Veltman 
scheme and the DRED scheme used by the Italian group in 1981. I am sure that presently a 
calculation of that type is done fully automatically by means of appropriate computer programs 
but in the winter 1988/1989 and in the spring of 1989 we did it almost entirely by hand. This 
has the advantage that each step of the calculation can be followed and enjoyed in contrast to 
looking constantly at the computer. 

Working then entirely by hand the first step is the calculation of two-loop momentum in- 
tegrals keeping 1/e^ and 1/e terms. In doing this, it is useful to factor out any Dirac and 
colour structures. In this manner the results for the two-loop integrals in question can be 
used for the calculations of two-loop anomalous dimensions of other operators and in different 
renormalization schemes in which 75 is treated differently. This first step is rather tedious but 
straightforward and it is not surprising that working independently we obtained the same results 
for all 28 diagrams already in the first comparison. I have used these results for the calculation 
of two-loop anomalous dimensions of other operators in 1991 and 1999. I will discuss these 
calculations later on. In calculating colour factors I found the paper j^89j very useful. 

The last step of the calculation, the manipulation of Dirac structures turned out to be the 
crucial part of our work. When calculating the subdiagrams in Fig. [2] we encountered structures 
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(a) (b) (c) 



Figure 2: One loop current-current diagrams that contribute to one- loop anomalous dimensions 
and enter as subdiagrams the two-loop calculations of two-loop anomalous dimensions. The 
4- vertex "(g) denotes the insertion of a 4-fermion operator Qi. 



like 

r.7p7M <^ ^"Yl^, = 7,(1 - 75) (3.5) 

that we had to reduce to the operators Qi and Q2 which have the structure F (8) F. At one- 
loop one can do this in D = 4 dimensions as 1/e is the leading singularity, but in a two- loop 
calculation the 1/e singularity, from which the anomalous dimensions are extracted, is next-to- 
leading. Consequently 0{e) terms in Dirac structures multiplied by the leading l/e"^ singularity 
from the momentum integrals have an impact on two-loop anomalous dimensions and have to 
be taken properly into account. 

As the first scheme we considered the one with anticommunicating 75 in D 7^ 4 dimensions 
giving it the name NDR (naive dimensional regularization) . In order to reduce the structures 
like the one in p.Sp to Qi and Q2 we first followed the procedure of Tracas and Vlachos [10] 
who simply wrote 

r.7p7M <S> T'^Yl'' = AT,(g>T'' (3.6) 

and found the coefficient A by replacing by 7r and contracting the Dirac indices on both 
sides. This procedure, to be called the "Greek method" in what follows, gives A = 4(4 — e) and 
consequently 

r.7p7M ^ ^"Yl^ = 4(4 - e)F, ® F^ (3.7) 

This method is very efficient and can be applied to Dirac structures with many 7^ that 
appear at two-loop level. It can also be easily generalized to other operators. For instance in 
the case of 7^(1 — 75) (8) 7'^(1 -|- 75) one should replace (8) by 1. 

Applying this method to 28 diagrams in question and using it also for the counter diagrams 
we could readily find the total 1/e singularity. Subsequently including the two- loop wave function 
renormalization for the external quarks we found the two-loop anomalous dimension matrix in 



28 



the basis 



Q+ 



Q2 + Q1 
2 



Q 



Q2-Q1 



2 



(3.8) 



We expected this matrix to be diagonal but we both found it to contain non-diagonal terms. 
There was a hope for an hour. We disagreed on four diagrams. After a fortunate draw 2 : 2 
we were certain that our calculation was algebraically correct but unfortunately the unwanted 
non-diagonal terms were still there. The only place, we could think, something went wrong 
was the "Greek method" for the reduction of the complicated Dirac structures to the physical 
operators Qi and Q2 as given in (j3.7p . 

Indeed in addition to the physical operators Qi^2 on the r.h.s of (13. 7p one could have other 
operators with different Dirac structures that vanish in D = 4. They had to vanish in D = 4 
because in D = 4 the formula (j3.7p is correct as can be easily checked by using standard 
manipulations of Dirac matrices. Peter called these operators "effervescent" operators. I did 
not object to this name. How could I? After all he is English, not me. Later we were told that 
the proper name is "evanescent" and we used this name in the following papers but in our first 
paper published in Nucl. Phys. B still "effervescent" operators appear [90]. Amusingly our 
Rome competitors used our wording still several years later [181]. 

Explicit expressions for the evanescent operators can be worked out. We have done it in our 
paper. But in practice it is more convenient to define them simply as the difference between 
the r.h.s. and l.h.s. of (|3.7p and to insert them like that in the relevant two-loop diagrams. 
Therefore instead of (|3.7p we have 



with the evanescent operator E^^^ defined simply by this equation. As discussed in detail in 
[9H |92] this is not the only possible definition of evanescent operators but possibly the most 
convenient one. The unphysical arbitrariness in the definition of evanescent operators has also 
been emphasized by Jamin and Pich [109J. 

Having indentified the possible origin of our problems we have incorporated the evanescent 
operators into our calculation. In particular we derived, to my knowledge for the first time, 
formulae that allow the extraction of the two-loop anomalous dimensions of physical operators 
in the presence of evanescent operators. The outcome of these efforts is section 4 of our paper. 
We have invested plenty of time in writing this section but apparently several of my colleagues 
had a difficult time in following it. I tried to improve on it in my Les Houches lectures [16] . 
where a systematic procedure for the inclusion of evanescent operators in the calculation of two 
loop anomalous dimensions of local operators can be found. 

Having the full machinery for the evaluation of the contributions of evanescent operators at 
hand, we could now find that their presence not only modified the diagonal terms in the 2x2 
matrix in question but also canceled the off-diagonal terms. We were now sure that the first 



r7p7M ® ^7^7^ = 4(4 -e)r0r + E- 



iNDR 



(3.9) 
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(28) (29) (30) 



Figure 3: Two-loop current -current diagrams contributing to 7s . The curled lines denote 
gluons. The 4-vertices "(8> (8>" denote standard operator insertions. In addition shaded blobs 
stand for self-energy insertions. Possible left-right or up-down reflected diagrams are not shown. 

NLO summit has been conquered: we knew the or equivalently the {Qi,Q2) matrix 

at the two-loop level in the NDR scheme. As at no place in the calculation it was necessary 
to evaluate the dangerous traces Tr^j^'jp'y^'yx'y^), we were quite confident that this result was 
correct. 

The calculation in the 't Hooft-Veltman (HV) scheme for 75 was technically more difficult 
because of the horrible Dirac algebra for which we had to use a computer program written by 
Peter. As I did the NDR calculation entirely by hand, including Dirac algebra, I could test the 
correctness of this program. Otherwise we did not encounter any obstacles and we soon had 
the two-loop anomalous dimension matrix of {Qi,Q2) in the HV scheme. By calculating the 
one- loop diagrams in Fig. 2 we found the matrix Af in ()2.40p . relating the matrix elements of the 
operators in question in the NDR and HV schemes. Inserting it in (|2.39|) we could indeed verify 
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that our results for the two-loop anomalous dimensions in these two schemes were compatible 
with each other. 

Strictly speaking this was the end of the story as we had the two-loop anomalous dimension 
matrix in the HV scheme that did not have any mathematical inconsistencies related to 75. We 
have also demonstrated that at least in this case a consistent calculation in the simpler NDR 
scheme could be made. Still we were curious whether the calculation of the Italian pioneers in 
1981 was compatible with our results. Actually it would suffice to calculate the relevant matrix 
Af relating HV or NDR scheme to the DRED scheme to find out that the calculation of 1981 
was correct, but for reasons that I do not understand today, we repeated the 1981 calculation 
confirming diagram by diagram the results of the Italian team. Most probably we were feeling 
very strong in such calculations and we were simply delighted in producing results for these 28 
diagrams in a different renormalization scheme. 

We have submitted our paper [SD] to Nucl. Phys. B in June 1989 and sent our preprint to 
CERN and SLAC libraries. One should recall that in 1989 the Los Alamos archive did not exist 
yet. Few weeks later our preprint has been distributed by ordinary mail around the world. 

At the end of July 1989 I attended the Photon-Lepton conference in Stanford. On the first 
day of the conference the two Guidos of the 1981 team congratulated me because of our paper. 
They were truly delighted. They were apparently not sure that their paper was correct. Indeed 
the calculation in the DRED scheme is very involved as also some aspects of QCD coupling 
renormalization have to be modified. 

At the same conference I met for the first time Matthias Jamin, who just got his PhD in 
Heidelberg and was supposed to join my group in Munich two months later. I told him about 
the MNLC and he immediately became the third member of the club. Already in October of the 
same year we were climbing together the second NLO summit, the QCD corrections to AF = 2 
processes that I will describe briefly in the next section. For the time being I will continue 
with the AF = 1 effective Hamiltonian for non-leptonic decays including now the QCD and 
electroweak penguin operators. 

In Table [2] we collect the references to the papers which calculated NLO and NNLO cor- 
rections to AF = 1 processes except for rare and radiative decays discussed in Sections 6-8. 
Two-Body B Decays in QCD Factorization (QCDF) approach are discussed in Section 9. I thank 
Gerhard Buchalla for helping me in collecting the references to NLO and NNLO calculations in 
QCDF given in this table. 

3.3 QCD Penguin Operators 

In the fall of 1990, after two successful expeditions, it was time to return to the QCD penguin 
operators that were the main topic of the seminal supper with Guido Martinelli in the Ringberg 
castle two and a half years before. There were no signs coming from Rome that the Italian 
team was making any progress on penguins but I started worrying that they were far ahead of 
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Table 2: NLO and NNLO Calculations for Non-leptonic and Semi-Leptonic AF = 1 Transitions 



Decay 


NLO 


NNLO 


Current-Current \Q\^Ll2) 
QCD penguins (Q3, (^4, Q5, Qe) 
electroweak penguins (Q?, Qs; Qq; Qio) 

inclusive non-leptonic decays 
Two-Body B-Decays in QCDF 
Current-Current (BSM) 
Penguins (BSM) 

Semi-Leptonic B Decays ( Vcfe , Ktfe ) 


[51190] 

[Ml ESI ESI Ea EH EH], [100] 
[inilESlEaEl] 

|1U21 11031110411105] 

|5U6l 11061 11071 1108): |109j 

|279]-|289] 

[301 [29] 

[29] 

|121j-|126|.|134lll35| 


[93] 

m 

[S3] 

|295|-|303j 

|126|-|133j.|136j-|139j 



us. For this reason I decided to increase our team. Markus Lautenbacher, my former diploma 
student and since April 1990 my PhD student, became the fourth member of the MNLC and its 
first PhD student. Markus did not have any experience with two-loop calculations but his high 
computer skills and an impressive discipline in doing research convinced me that he would be 
a great help in our project. Our first goal was the calculation of the 6x6 two- loop anomalous 
dimension matrix 7^^^ describing the mixing under renormalization of the operators Qi, Q2-, --Qe 
in the NDR and HV schemes. The calculation of ji^^ involves the insertions of all these operators 
into vertex diagrams considered already by Peter and myself in our first paper and into two-loop 
penguin diagrams in Fig. |3]to be considered for the first time. The latter diagrams do not have 
any impact on the sector {Qi, Q2) so that the corresponding 2x2 submatrix of 7^^^ calculated 
by Peter and myself remained untouched. 

In the first month I worked closely with Markus helping him in making first steps on this 
new ground, whereas Peter and Matthias worked independently by themselves. Later Matthias 
and Markus worked closely together and constructed an efficient program for Dirac algebra ma- 
nipulations in 7^ 4 in the NDR and HV schemes [110] . This program written in Mathematica 
became an important part of our project in particular in the case of the HV scheme and even I 
used it in spite of my previous comments on computer manipulations. In this scheme the calcu- 
lations of two-loop penguin diagrams by hand were prohibitive as even computer manipulations 
required in 1990 a good PC. The evaluation of two-loop momentum integrals in a few penguin 
diagrams that we did mostly by hand turned out to be rather involved and a method by Peter to 
find the coefficients of the divergences in these particular diagrams numerically was very helpful. 
The corresponding calculation of the vertex diagrams was simple as we had already all integrals 
from |90] . Only the Dirac structures were different. 

There were two new features with respect to the calculation of the {Qi,Q2) system of 1989. 
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First we had to face the dangerous traces Tr (7^ 7p 7,^7^75). In the HV scheme they can be 
straighforwardly evaluated but their evaluation in the NDR scheme could lead to wrong results. 
For a few weeks we thought that we had to abandon the calculation in the NDR scheme but at the 
end we solved the problem in two ways. My solution was to work, dependently on the diagram 
considered, with a second operator basis {Qi} with Qi being Fierz transformed operators of Qi. 
With the help of these operators it was possible to avoid the appearance of the dangerous traces. 
However, simply replacing Qi by Qi in order to avoid dangerous traces with 75 and inserting 
it into a two-loop penguin diagram would eventually give the wrong result for 7^^''. This is the 
second new feature of the presence of penguin diagrams: the insertion of (Qi, (32, Qs, Q4) into a 
two-loop diagram gives different result from the insertion of {Qi,Q2,Q3,Q4) and consequently 
the two-loop mixing between {Qi,Q2) and {Q3,Q4) differs from the one between {Qi,Q2) and 
{QsjQi). Fortunately, similary to the renormalization scheme dependence of the difference 
in question could be found by performing a one-loop calculation that did not involve dangerous 
traces with 75. Incorporating this difference properly into the calculation that involved both the 
original operators and their Fierz transforms allowed then to obtain 7^^"* for the original basis 
in the NDR scheme without any problems with 75. 

This procedure is described in detail in |95j . In the case of the evaluation of the two-loop 
anomalous dimensions of electroweak penguin operators it had to be generalized to include more 
bases because of the more complicated flavour structure of these operators. This is described 
in |101j . This procedure was followed by Matthias and Markus. Peter succeeded to avoid the 
dangerous traces in a different manner but I do not remember how. This is however immaterial 
as our independent results for 7^^^ in the NDR scheme agreed with each other. 

I must admit that at that time I was very satisfied with my procedure of working simul- 
taneously in D ^ 4 with the original basis and the Fierz transformed basis and making finite 
renormalizations through one-loop calculations at the end as explained above. However, in 1994 
a more elegant and a more systematic procedure with the same outcome has been proposed by 
Matthias Jamin and Toni Pich |109j in the course of their NLO analysis of inclusive AF = 1 
transitions. Basically one can make four dimensional Fierz transformations in a 7^ 4 calcula- 
tion provided the evanescent operators that vanish in D = 4 under the Fierz transformation are 
also included in the analysis. They are simply given by the difference of a given operator and 
its Fierz transformed operator. This method has been rediscovered by Mikolaj Misiak and Jorg 
Urban in the process of another NLO climb that I did with them in 1999. I will return to it in 
Section 3.5. 

Still my cooking recipe that involved calculating the differences between the one-loop inser- 
tion of an operator and of its Fierz transformed operator is very useful for finding out whether 
the evanescent operators of that type are relevant for the calculation of two-loop anomalous 
dimensions or not. If this difference vanishes, the evanescent operators in question do not con- 
tribute. This is the case of all operators Qi,-.Qio inserted into the current-current diagrams 
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Figure 4: Two-loop penguin diagrams contributing to 7^ . The curled lines denote gluons. 
Square- vertices stand for penguin insertions. Possible left-right reflected diagrams are not shown. 

and of Q5, ..Qs inserted into the penguin diagrams. 

We thus had in the Spring of 1991 the full 6x6 anomalous dimension matrix 7^ at 0{a1) 
in the NDR scheme^ The calculation in the HV scheme was time consuming because of the 
difficult Dirac algebra but with the computer program developed by Markus and Matthias we 
could calculate 7^^^ in the HV scheme both directly by calculating the traces with 75 and by 
using my procedure discussed above, obtaining the same result. Finally calculating the relevant 
one loop shift Af we verified that our results for 7^^^ in the NDR and HV schemes were consistent 
with each other. 

Next we calculated the initial conditions for the Wilson coefficient functions at /i = 0{Mw) 
both in the NDR and the HV scheme and verified that the scheme dependence of these coefficients 

'it's generalization to 10 x 10 matrix is discussed in the next subsection. 
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cancelled the one of the evolution matrices at /x = 0{Mw) as explained in section 2. 

Thus the third NLO summit has been reached. As usually reaching a summit a photo is 
taken, I thought it was appropriate to show the resulting two-loop anomalous dimension matrix 
7s^^ in the NDR scheme: 



/ _2i _ 2/ 

/ 2 9 



7 + 2/ 
2^3 









2/ 
3 



21 _ 2/ 
2 9 









79 
9 

202 
243 



5911 
486 



+ 



71/ 



_7 
3 

1354 
81 



5983 _|_ / 
162 3 



_65 
9 

1192 
243 



_7 
3 

904 
81 



379 I 56/ 
18 243 



-61/ 
9 



.91 + 



808/ 
81 



2384 
' 243 

130 _ 
9 



v 











-682/ 
243 



-11/ 



106/ 
81 



71 



+ 



225 I 
2 



_ 71/ 
9 

502/ 
243 

61/ 



1676/ 
243 



1808 
81 



14 I 646/ 
3 81 



-99 + 



11/ 



(3.10) 



1343 _|_ 1348/ I 
6 81 / 



This matrix looks truly horrible and I will in what follows refrain from showing other photos of 
this type. In fact after the inclusion of electroweak penguins and going to three loops the results, 
although very impressive, cannot be easily digested. On the other hand the 2x2 submatrix in 
the upper left corner, our first summit, conquered by Peter Weisz and myself in June of 1989 
looks beautiful and simple. 

Our paper has been submitted to Nucl. Phys. B in May 1991 [94]. In addition to the 
6x6 matrices and the Wilson coefficients of {Qi,Q2, ■■■Qe) in the NDR and HV schemes for 
both AS" = 1 and AB = 1 decays it contained general expressions for the evolution matrices 
U{fii,fi2) including NLO corrections. In order to derive these expressions we have used the 
general all order formulae of my 1980 review on asymptotic freedom in deep inelastic scattering 
[lllj . Finally we have demonstrated the scheme independence of the resulting decay amplitudes. 
Thus at last, three years after the Ringberg workshop, the Wilson coefficients of current-current 
and QCD penguin operators were known at NLO in the NDR and HV schemes. 

I have presented these results in a parallel session at the joined Photon-Lepton and European 
Physical Society Meeting that in 1991 took place in Geneva, Switzerland |112j . Rather disap- 
pointigly only few of my colleagues appreciated these results. In particular Eduardo de Rafael 
thought it was an overkill in view of the uncertainties in the hadronic matrix elements of the 
operators in question. Eduardo got interested in our work only ten years later, when he wanted 
to know more about the role of evanescent operators and Fierz relations in our calculations that 
he and his collaborators wanted to combine with their calculations of hadronic matrix elements 
within the large N approach. 

Also to my great surprise there was essentially no reaction from Guido Martinelli. He only 
informed me that his PhD students are working on this project as well and that in order to 
complete the project one needs the hadronic matrix elements of QCD penguin operators. In 
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fact in 1991 I knew these matrix elements in the context of 1/A^ expansion developed with 
Bardeen and Gerard several years earlier but Guido meant here the ones obtained by lattice 
methods. Unfortunately Guido's dream to achieve the latter result did not materialize until 
today. Hopefully in this decade we will know whether lattice results for the QCD-penguin 
hadronic parameter Bq will be very close to large-N values as it turned out to be the case of Bk 
relevant for the CP-violating parameter ex in — ?• vrvr {113j. 

As far as QCD penguins are concerned the first NLO analysis of penguin induced i?-decays 
and the related CP-asymmetries using our two-loop results has been performed by my student, 
Robert Fleischer, in the Summer of 1992 |lUUj . Robert combined his one-loop calculations of 
matrix elements with the two-loop anomalous dimensions discussed above and demonstrated 
the scheme independence of the final result. It was his Diploma thesis. During his Phd studies 
Robert fall in love with electroweak penguins and investigated their role in non-leptonic decays. 
Therefore he did not have time to participate in the subsequent papers on NLO QCD corrections 
to weak decays except for his second paper in |100j which is a proof of his new interests in 1994. 

Further progress in evaluation of the 6x6 anomalous dimension matrix for the operators 
Qi, . . . , Qq will be described in Sees. 16.3.31 and 16.41 

3.4 Electroweak Penguin Operators 

In spite of this rather moderate interest in our work in 1991 I was convinced that we should 
continue our project. The next step was to extend our calculation of ji^^ in (|2.37p to electroweak 
penguin operators Qj, Qs, Qg and Qio and to calculate the ten dimensional two-loop anomalous 
dimension matrix ^se^ that is necessary for the inclusion of the electroweak penguin operators 
at the NLO level with the goal to calculate the CP- violating ratio e' /e. Moreover, we wanted 
to write up the details of all these calculations. We have not done this in |94j . 

Unfortunately, Peter told me that he would only be involved in the calculation of 7^^^ as he 
was again very much involved in the collaboration with Martin Liischer. Thus only Matthias, 
Markus and me were involved in the 7^6'' project that amounted in particular to the calculation 
of the two-loop diagrams in Fig. [5l Having all the machinery at hand we performed both 
calculations during the fall of 1991 and the winter 1991/1992 so that in March 1992 we had 7^^^ 
and jil^ including all the ten operators in the NDR and HV schemes. Moreover we calculated 
0{a) corrections to the Wilson coefficients at fj, = My/, an ingredient of the NLO analysis, that 
is necessary to remove the renormalization scheme dependence from the decay amplitudes. 

Unfortunately, there was a problem with our result for T-iP. While the Yi^P]ndr and [75^^]^!/ 
were compatible with each other, Yi^)]ndr and [7es^]Hy were not. That is [t'IP]//!/ obtained by 
the direct two-loop calculation differed by a small amount from [7es^]_ffv found from Yi^)]ndr 
by means of a formula analogous to (j2.39p that is given in |101j . 

We spent some time in order to clarify this discrepancy but after a few weeks we made a 
pause in our search for the error. I think we were simply exhausted. Moreover everyone was 
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Figure 5: Two-loop current -current diagrams contributing to ^se . The wavy lines denote gluons 
or photons. The 4-vertices "(8> (8>" denote standard operator insertions. Possible left-right or 
up-down reflected diagrams are not shown. 

involved simultaneously in other projects: Peter with Martin Liischer, Matthias and Markus 
with writing up their paper on "TRACER", the program, written in Mathematica, for Dirac 
algebra va. D ^ A for NDR and HV schemes |110j and in my case in addition to being the head 
of the theory institute at our university and finishing a review article with Michaela Harlander, 
I started a new NLO climb: in the spring of 1992 our NLO club got the fifth member, a very 
important one, namely my new PhD student Gerhard Buchalla with whom I planned to attack 
at the NLO level all rare semi-leptonic K and B decays dominated by Z'^-penguins. More about 
this in section 5. 

Fortunately before making a pause in our climb we decided to write up the two papers, one 
including Peter on 7^^^ that did not have any problems and the second one without him on 
that had the problem mentioned above. Thus already in April 1992 our papers were essentially 
finished but before we could show them to the public we had still to solve the remaining problem 
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in the second paper. 

From the beginning I was fully confident that our results in the NDR scheme were correct. 
The calculations were simpler than in the HV scheme and I was able to make several tests that 
all worked. As the calculations of Ar matrices relating the NDR and HV schemes is a one loop 
affair I was also confident that our results in the HV scheme obtained from the NDR scheme by 
means of a relation similar to (j2.39p were also correct. However, the game of making a given 
NLO calculation in various schemes and checking the compatibility of the results started by Peter 
and myself three years before somehow fooled us and we did not send the papers for publication 
although we had all results already in April 1992. In the language of moutain climbing it is 
afterall irrelevant whether the first climb of a summit was done using the NDR "climbing style" 
or HV one. 

Fortunately, the Rome group did not present their results at the 1992 summer conferences 
and consequently we were still in the game. Moreover, Matthias became the CERN fellow and 
could inform us in the first days of November 1992 that Guido Martinelli will give a seminar 
on e' je beyond leading logarithms four weeks later. It was time to be active again. Feeling 
like colonel Hunt before the final attack to conquer the Mount Everest summit I convinced my 
collaborators to send out the paper on 7!^^ in the existing form and to present the details of the 
calculation of 7es^ in the second paper only in the NDR scheme making the shift Ar to obtain 
it in the HV scheme. Our two papers |951 llOlj appeared in the second half of November 1992, 
roughly two weeks before Guido's CERN talk and three weeks before the Rome group sent out 
their letter to the Los Alamos archive |97] . 

To our delight the Rome team consisting of Marco Ciuchini, Enrico Franco, Guido Martinelli 
and Laura Reina agreed with our results on the anomalous dimension matrices in both NDR 
and HV schemes but to our surprise they did not present any details of their calculations of 
these matrices. Instead they presented their analysis of e' je including NLO QCD and QED 
corrections. Thus although the Munich team has published as the first group all ingredients of a 
NLO analysis of IS.F = 1 processes: two-loop anomalous dimensions and the Wilson coefficients 
at /i = Miy, the Rome group was the first to present a NLO analysis of e'/e that included both 
QCD and electroweak penguin contributions. 

Our NLO analysis of e' je appeared in March 1993 in a very long paper (112 pages) that 
fortunately was accepted in this form in Nucl. Phys. B. [96]. We have presented there very 
explicit formulae for the Wilson coefficients of all operators and studied also the /i dependence 
of the hadronic Bi parameters. Moreover, we have proposed a method for extracting some of 
these parameters from the CP-conserving amplitudes. In this manner we could incorporate the 
A/ = 1/2 rule into the analysis of e' je. The results of this analysis have been used by us in the 
1990's to gradually develop an approximate but rather accurate formula for e' je that depends 
on three parameters iZg) -^8 and O/^ with the first two representing the relevant matrix element 
of and Qs operators and the last one summarizing isospin breaking corrections [1171 1118] , 
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Figure 6: Two-loop penguin diagrams contributing to jse . The wavy lines (denote gluons or 
photons with each diagram containing one gluon and one photon. Square- vertices stand for two 
types of penguin insertions. Possible left-right reflected diagrams are not shown. 

The last update of this formula can be found in |119| . 

As the formulae presented by the Rome group were not as explicit as ours, an analytic 
comparison between their final results and ours was not possible but numerically the agreement 
was within a few percent. 

In our long paper ^96j on e'/e we have also presented for the first time NLO Wilson coefficients 
of the operators relevant for non-leptonic AB = 1 decays. Amusingly the interest of particle 
physics community in e' /e, in the middle of 1990's, was rather moderate and our e'/e paper of 
1993 was cited in this period mainly for our B physics results. 

One month after our final paper of this period the Rome group presented in great details 
their calculation of two-loop anomalous dimensions concentrating on the HV scheme |98j. This 
helped us to identify the error in our direct calculation of [7es^]_f/y- Their result for this matrix 
indeed agreed with ours obtained in our paper indirectly through [jil^jNDR- 

The most important result of this Munich-Rome competition that lasted for several years is 
the agreement on the two- loop anomalous dimensions of the operators Qi{i = 1, ...10) that have 
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been used by the flavour community since then. While they were calculated with the aim to 
find e'/e at the NLO level, they play also essential role in all AF = 1 transitions, not only the 
non-leptonic ones. Thus they also enter the NLO calculations of B ^ Xg'j, B — t- X^e^e" and 
even K — )• ttvU decays. This is why in planning the grand expedition in 1988 it was essential to 
calculate these anomalous dimensions first. 

In March 1993 I gave a seminar on our e'/e analysis in the theory group at CERN. I was 
approached there by two young physicists whom I did not meet before. It was Marco Ciuchini 
and Enrico Franco who came from Rome on a night train for a day to CERN to listen to my 
talk. I was really impressed that our competition went so far. However, it turned out that this 
friendly competition did not end in 1993. 

3.5 More Operators 

In 1997 an Italian group consisting of six climbers [30], including three from the expedition of 
early 1990's calculated two-loop anomalous dimensions of a set of operators relevant at NLO 
for AVg ii in — B^^i mixings and in particular for AF = 1 and AF = 2 non-leptonic decays 
and transitions in the extensions of the SM, like MSSM at large tan/3, multi-Higgs models and 
generally models that include in addition to left-handed currents also right-handed currents and 
scalar currents. These operators are given in the operator basis of [25] in (j2.15p and (j2.16p . 

I was a bit surprised that we have not looked at these operators, except for the first one, in 
Munich before. Otherwise we would calculate their anomalous dimensions already in 1993. The 
whole machinery developed by us at the beginning of the 1990's could be used here. Knowing 
this, there was essentially no point in repeating this new Italian calculation. The results given 
in [30] were bound to be correct. However, one result presented in this paper and even pointed 
out in the abstract made me interested in looking at this analysis closer. 

The Italian group calculated the two-loop anomalous dimensions of the operators in ()2.15p 
and (|2.16p in the so-called RI scheme that is apparently useful for lattice calculations and in 
the NDR scheme. The calculation in the RI scheme was entirely new. The calculation for the 
operators Q12 in the NDR scheme was really not new as the anomalous dimensions of these 
operators can be directly obtained from our earlier calculations of QCD penguin operators. 
What was new in the NDR scheme were the two-loop anomalous dimensions of the operators 
Qi2^. Here the Italian group found a surprising result: the analogs of Q+ and Q_ operators in 
()3.8p mixed in this sector under renormalization. Their two-loop anomalous dimension matrix 
was non-diagonal and the results for the NDR presented in the appendix of this paper looked 
very complicated. If this was indeed true, the NDR scheme as defined by Peter and myself in 
|90| would not be an elegant scheme. 

Also Mikolaj Misiak, who joined MNLC in 1995, was interested in this result. Together 
with Jorg Urban, a PostDoc in my group, we decided to look at the Italian paper closer. I 
should have noticed it right away, but it was Mikolaj who reminded me of the "Greek story" 
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of missing evanescent operators in the case of Qi and Q2 and the resulting mixing of and 

found by Peter and myself in 1989. This time the Fierz vanishing evanescent operators 
were involved. I have discussed them already in section 3.1. We soon suspected that Italian 
masters were performing Fierz transformations in a D 7^ 4 calculation without including the 
evanescent operators in question. We have all three independently performed the calculation of 
the anomalous dimensions of this sector now including the Fierz vanishing evanescent operators 
obtaining the result in the NDR scheme that was much simpler than the one of the Italian 
group, in particular the non-diagonal entries disappeared as they did in 1989 in the case of the 
current-current operators (Q+jQ-)- 

Mikolaj and Jorg confirmed the RI calculation of [30] in an arbitrary covariant gauge and 
found the matrices Ar relating the RI scheme and the NDR scheme. In this manner we could 
also find the compatibility of the RI result in [30] and our NDR result but the one loop matrices 
connecting these two schemes were clearly different from those given in |30| . 

Unfortunately, our Italian colleagues did not agree with our interpretation of the strange 
form of their result but they admitted that the renormalization scheme they used was not the 
standard NDR scheme of [90] and that our result was more elegant and phenomenologically 
more useful. In summary: the first NLO climb in the RI scheme related to the operators (|2.15|) 
and ()2.16p should be credited to the Italian group while in the NDR scheme to us. 

The operators in (j2.15p and (j2.16p do not constitute the full set of six-dimensional four quark 
operators contributing to = 1 processes. In addition to QCD penguins and electroweak 
penguins of the SM there are other penguin operators. In our paper [29] we have therefore 
generalized our analysis of two-loop anomalous dimensions to the full set of AF = 1 four-quark 
operators. These results are much less known but should be useful in the extensions of the SM 
one day. 

3.6 QCD Corrections to Semi-leptonic B Decays 

This is probably a good place to summarize the QCD corrections to semi-leptonic B decays that 
are necessary for an accurate determination of the CKM elements \ Vcb\ and \ Vub\- As I did not 
take part in these calculations I consulted a very prominent member of the MNLC club, Paolo 
Gambino, who contributed in an important manner to these calculations. As we will soon see 
another member of our club, Andrzej Czarnecki, who will enter the scene later on, also made 
important contributions to this field. A nice summary of this topic as far as 6 — ?■ civ is concerned 
can be found in [120] . 

In what follows I will list papers where only the purely perturbative corrections to semilep- 
tonic B decays have been computed leaving the summary of power corrections to other reviewers. 
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3.6.1 NLO Corrections to Inclusive B — )• XJi^ Decays 

The pioneering first steps in tliis field can be found in |121j . Tlie first complete analytic calcula- 
tions, not only rate but also a few differential distributions, have been performed much later by 
Andrzej Czarnecki, Marek Jezabek and Hans Kiihn |122) . In this context also the corrections to 
the rate in a compact form found by Yossi Nir should be mentioned |123j . 

Next corrections to moments of hadronic spectra which used the results of previous calcula- 
tions can be found in |124| . while full triple differential distribution at 0{as) that are necessary 
for realistic experiments have been obtained in [125] and in particular in [126J. 

3.6.2 NNLO Corrections to Inclusive B — > XJf Decays 

First BLM-NNLO corrections to the rate have been obtained in |127j and in particular in |128j . 
The BLM-C'(q;^/3o) corrections to lepton spectrum can be found in [129] and to to triple differ- 
ential distributions in |126j . 

Next non-BLM two-loop corrections (analytic, zero cut) have been obtained in |130j and 
in the numerical form but with realistic cuts in [131j . Finally complete two-loop corrections 
at specific kinematic points (zero recoil, which is important for the determination of \Vcb\ from 
B D^*hiy) can be found in [l32]. 

3.6.3 NLO and NNLO Corrections to Inclusive B XJui Decays 

NNLO complete calculation of the width has been performed in |133j . NLO and BLM-NLO 
full triple differential distributions have been calculated in |134j and |135j . respectively. Leading 
shape functions and resummation in B — )• Xulfi has been done in |136j and non-leading shape 
functions in |137j . The most recent calculations dealing with this topic are |1381 1139j . 



4 AS' = 2, AB = 2 and AB = Transitions 

4.1 Effective Hamiltonians for Ai^ = 2 Transitions 

Let us begin this section by recalling the effective Hamiltonians for A5 = 2 and Ai? = 2 
transitions in the SM. We have first [HDl [HH 



a/A5=2 
^eff 
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Q{AS = 2) + h.c. 



where A, = V*gV-^. Here /i < /ic = 0{mc). In (j4.ip . the relevant operator 
Q{AS = 2) = {sd)v-A{sd)v-A 



(4.1) 



(4.2) 
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is multiplied by the corresponding coefficient function. This function is decomposed into a 
charm-, a top- and a mixed charm-top contribution with So{xi) and So{xi,Xj) being one-loop 
box functions in the SM. 

Short-distance QCD effects are described through the correction factors rji, r]2, Vs the 
explicitly ag-dependent terms in the last line of (j4.1|) . This factor allows to introduce the 
renormalization group invariant parameter Bk by 



-2/9 



(4.3) 



(K°|(s(i)y-A(sd)y_A|i^°) = -BMP^ml. 



(4.4) 



The corresponding Hamiltonian for B^ ^ — B^ ^ mixing has similar structure but it is simpler 



as only the top contribution matters. We have for B^ — B^ mixing 



njAB=2 
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Here Hb = 0{mb), 

Qi{AB = 2) = {bq)v^Aik)v-A, q = d,s. 



Q'^iAB = 2) + h.c. 



(4.5) 



(4.6) 



The renormalization group invariant parameters Bq are defined by 
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(4.7) 



{B^q\{bq)v-A{bq)v-A\B^g) ^ -BB,{lj)Flm\ , 

where -Fb, is the i?g-meson decay constant. 

We are now ready to discuss the history of the NLO QCD calculations of 



(4.8) 



Vi ^ Vcc, 



(4.9) 



4.2 The Top Quark Contributions 

In the fall of 1989 Matthias Jamin joined my group becoming the third member of the MNLC. 
Instead of continuing our NLO calculations for AF = 1 transitions we decided to calculate 
the NLO QCD corrections to top quark contributions to the effective Hamiltonians for B'^^ — 
B^g mixings and — mixing. We were not the first to do this climb but the first two 
attempts were unsuccessful. The calculations were plainly wrong with the results for the Wilson 
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coefficients exhibiting infrared regulator dependence that was a consequence of an incorrect 
matching of the full and effective theories. Moreover these calculations did not include the two- 
loop anomalous dimension of the operator Q{^F = 2). As I have a high respect for the leaders 
of these two expeditions, that otherwise had significant contributions to our field, I prefer not 
to refer to these papers. 
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Figure 7: Examples of two-loop diagrams contributing to B^^ — mixing. 

Our new NLO project |142j involved two parts. One very easy, the other one rather difficult. 
The first one was the calculation of the two-loop anomalous dimension of the relevant operator 
Q{AF = 2). It took no time as the anomalous dimension of Q{AF = 2) is simply equal to 7+, 
the anomalous dimension of (5+ in (j3.8p calculated at the two-loop level by Peter and myself 
half a year before. The second part was much more time consuming. It involved the calculation 
of 0{as) QCD corrections to the box diagrams involving internal bosons, the Goldstone 
bosons ip^ and the top quark exchanges. This is a two-loop calculation in a full theory with 
massive heavy particles. Two examples of contributing diagrams are shown in Fig. [71 The 
remaining diagrams can be found in Fig. 2 of |142j . There are 8 classes of diagrams in total 
shown in this figure from which the remaining diagrams can easily be obtained. 

In order to extract the relevant Wilson coefficient of the operator Q(AF = 2) from this 
calculation the proper matching to the effective theory has to be made. This requires the 
calculation of the 0{as) corrections to the matrix element of this operator between the external 
quark states. The two-loop calculation of the 0{as) corrections to the box diagrams involves 
infrared divergences. We decided to set the external momenta to zero and regulate the infrared 
divergences by the masses of the external quarks. Working off-shell introduces necessarily gauge 
dependence (gluon propagator) in the final result. The final result of the box diagram calculation 
was therefore gauge and infrared regulator dependent. These dependences certainly did not 
belong to the Wilson coefficient of the operator Q{AF = 2) but to its matrix element and were 
removed in the process of matching of the full theory to the effective theory with the latter 
exhibiting precisely the same gauge and infrared regulator dependences. 

The main results of our paper were the values of the QCD factors rjB and r/2 for B*^^ — B^^ 
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and K — K Hamiltonians, respectively: 

T/B = 0.55 ±0.01, T/2 = 0.57 ±0.01 . (4.10) 

Our result for r^s has been confirmed several years later in |64j . 

These QCD factors have been used in the last 20 years frequently in the literature. While 
using these particular values one should remember that they should be used in conjunction with 
the renormalization group invariant parameters Bk and Bg^d in (|4.3p and (j4.7p . respectively. 
These parameters have been introduced at the LO already in 1983 by Wojtek Slominski, Herbert 
Steger, and myself |140j and have been generalized to NLO in the present calculation. I do not 
remember who is the father of the Bk parameter without the "hat", but I think it is John 
Donoghue. 

This takes care of the /x dependence at the lower end of the RG evolution. r]i resulting from 
the calculation of C(/x) and Bi representing the matrix element are so defined that they 

separately do not depend on /x. However, rji depend on the scale = 0{Mw) at which the 
matching between the full and the effective theory is made. This dependence is cancelled by the 
/iVK dependence of the initial conditions C{iiw) and this cancellation is only meaningfull at the 
NLO. More importantly the r]i are so defined that they multiply the leading order box functions 
and consequently they include the 0{as) corrections to the box diagrams that Peter, Matthias 
and me calculated. Consequently they depend also on the scale fit at which the running top 
quark mass mt{fJ,t) used in the calculation is defined. This fit dependence of r]i is important as 
it cancels up to higher order corrections the fit dependence of the Inami-Lim function So{xt{fit)) 
that remained uncompensated at LO. The values in ()4.10p correspond to fit = rrit. This turns 
out to be a convenient choice as with fit = mt, the QCD factors are practically independent 
of the actual measured value of mt{mt). The leading logarithm multiplying a large anomalous 
dimension of the mass operator vanishes at this scale. I mention this issue again because it is 
important. 

To my knowledge the issue of the fit uncertainty in /S.F = 2 transitions in the LO and its 
reduction through the inclusion of NLO QCD corrections, was for the first time addressed in 
|142j . I have described it here because it enters all the calculations presented below and our 
paper completed in the Spring of 1990 can be considered as the prototype of analogous two-loop 
calculations of the Wilson coefficients for any FCNC process that is sensitive to the top quark 
mass or any other heavy particle with colour, for instance squarks in the MSSM. 

I should emphasize that our calculation of NLO corrections to box diagrams with top quark 
exchanges performed in 1989-1990 had much simpler structure than the LO calculation of renor- 
malization group effects done by Fred Oilman and Mark Wise already in 1983 HB]- The reason 
is that in 1983 the typical values of mt considered in the literature were substantially lower 
than My/ and Fred and Mark in their 1983 LO calculation had to integrate out first and 
subsequently the top quark at much lower scales. In the range mt < fi < My/ they had to deal 
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with QCD corrections to bilocal operators originating in the contraction of the W propagator 
in a box diagram to a point but leaving the top quark propagator as it is. The renormahzation 
and the calculation of QCD corrections to these bilocal structures are rather involved even in 
the LO. Fortunately in 1989 it was already known that Myy ^ mt < 200 GeV. Consequently 
and the top quark could be integrated out simultaneously generating a local operator from 
the beginning and the QCD renormahzation of the bilocal structures at NLO, a formidable task, 
could be avoided in our calculation. 

Table 3: NLO and NNLO Calculations for AF = 2 and AF = Transitions 



Decay 


NLO 


NNLO 


m 
m 

ADMs BSM 


m 

|1421 [M! 
[I44lll45| 
[301 [29] 


[l46j 


AF = 2 Tree-Level 


[1501 11511 11521 1154| 
|152i I153i I154| 

m\ 





4.3 Charm and Top-Charm Contributions 

In the case of box diagrams with two charm quark exchanges there is no way out. One has to 
face the bilocal structures at NLO because the simultaneous integration of the charm quark and 
would lead to In Mw / n-c terms and consequently to the breakdown of perturbation theory. 
This rather difficult project was assigned in 1992 to my PhD student Stefan Herrlich. Stefan 
was a very good student and possibly he would succeed this climb by himself but fortunately 
for him and for the project, he was joined early 1993 by Ulrich Nierste. Ulrich got his diploma 
in Wiirzburg working with the loop masters like Manfred Bohm and Ansgar Denner and con- 
sequently he was fit for this difficult climb almost immediately after his arrival in Munich. In 
fact Ulrich was soon leading this important climb and developed to one of the most prominent 
members of the MNLC. 

The calculations of rji = rjcc and 773 = rjd QCD factors at NLO were completed in 1993 [143] 
and 1995 [144j . respectively and the full analysis of AS* = 2 Hamiltonian at the NLO level could 
be performed soon after [145j . These were truly heroic climbs that were not repeated by anybody 
until 2010, when two younger members of the MNLC performed the NNLO calculation of rj^: 
one of my many physics sons, Martin Gorbahn, a big star these days in multiloop calculations 
and my physics grandson, a PhD student of Uli Nierste, Joachim Brod [146j . It looks like the 
rji factors still remain in the possession of our physics family. 
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Even Guido Martinelli was impressed by the Herrlich-Nierste calculations and he told me 
this at least three times at different occasions. The values for i]i and r/3 enter the analysis of the 
CP-violating parameter ex and are relevant ingredients of any analysis of the unitarity triangle, 
in particular after the lattice value for Bk |113j and the estimate of long distance effects in 
ek improved recently [148t 1141^ . I was through all these years convinced, knowing Stefan and 
Ulrich, that the values for 771 and r/3 found by them were correct and the recent calculation of 773 
by Brod and Gorbahn confirmed my expectations. It would be good to know the NNLO value of 
r/i as well. Even with the technology of the present decade these NNLO calculations are heroic 
efforts. 

In fact Brod and Gorbahn completed this calculation few months after V2 of this review 
appeared |147j . The result is a bit disappointing as the NNLO calculation did not reduce the 
uncertainties in the charm part. Possibly with a better matching to the lattice calculations this 
can be achieved one day. 

4.4 Ar^ and ATd at NLO 

In the winter 1997/1998 I have been asked by Gerhard Buchalla, whether I would be interested in 
joining him, Martin Beneke, Christoph Greub, Alexander Lenz and Uli Nierste in the calculation 
of NLO QCD corrections to the life-time difference or equivalently ATg in the B'^ — B'^ system. 
At first I found it an interesting idea. Afterall ATg is much larger than AF^^ and working with 
my physics sons Alex and Uli, my physics stepson Martin and the Swiss master Christoph for 
the first time would be a real fun. In 1984 I studied AF^ at LO with Slominski and Steger jl4Tj 
and in addition after nine years of NLO climbing I was well prepared for this new expedition. 
Yet, this winter I was busy with writing up my Les Houches lectures and other project and I 
did not join Gerhard et al. Two papers in 1998 and 2002 |150^ 11531 1154) as well as numerous 
phenomenological analyses of Lenz and Nierste in the last decade resulted from this project 
|155| . Equally important this topic was the subject of the PhD thesis of Alexander Lenz. 

Interestingly, also the younger generation of Rome masters got involved in these efforts 
|15H 1152] . in particular Cecilia Tarantino, who became in 2006 my close one-loop collaborator 
within the Littlest Higgs Model with T-parity. I hope one day somebody will report on this 
competition, although as most of the authors are at least twenty years younger than me, it will 
still take some time. The results of these four papers played already an important role in the 
analyses of the Tevatron data and will certainly be very important for the LHCb. 

4.5 NLO QCD Corrections to Tree-Level AF = 2 Processes 

NLO QCD corrections to box diagrams are rather involved and it appears a bit premature to 
calculate them for extensions of the SM. Therefore, between 2006 and 2011, when I exclusively 
studied extensions of the SM, I decoupled from calculations of NLO QCD corrections. However, 
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in the fall of 2011 I noticed that for AF = 2 processes mediated at tree-level by colourless 
neutral gauge boson and neutral Higgs exchanges the matching conditions at NLO require only 
one-loop calculations and can be done model independently without too much effort. Somehow 
QCD experts including myself did not notice it before. 

In April 2011 Jennifer Girrbach, PhD student of Ulrich Nierste and consequently my grand- 
daughter in physics, joined my group. Jennifer had no experience in QCD calculations but I 
thought it would be fun to perform this climb with her and to teach her this field. Moreover, 
she was one of the stars of Ulrich's group and I was sure that we will reach this summit to- 
gether. Indeed Jennifer learned the QCD technology in a short time and performed all necessary 
calculations independently of me. We published our results on AF = 2 transitions already in 
January 2012 [156j . Since then we could use these results for concrete models with Z' tree-level 
exchanges. We could subsequently extend this calculation to non-leptonic AF = 1 processes 
mediated by colourless neutral gauge boson and neutral Higgs exchanges [157j . 



5 Rare K and B Decays 

5.1 Effective Hamiltonians for ir^uu and K^, — >■ tt^vv 

These decays have been with us already for more than twenty years and the world of particle 
physics is waiting for their precise measurements. A review published in 2008 |158j summarizes 
the status of two years ago. More recent developments are presented below and in |84j. 
The effective Hamiltonian for — )• tt^vv can be written as 

-H,^ = ^ " E (v:sVcdX\xc) + V:sVtdX{xt)){-sd)v-A{nyi)v-A. (5.1) 
V 2 27r sm^ t^w , ^ ^ 

The index l=e, ^, r denotes the lepton flavour. The dependence on the charged lepton mass 
resulting from the box-graph is negligible for the top contribution. In the charm sector this is 
the case only for the electron and the muon but not for the r-lepton. 
The function X{xt) relevant for the top part is given by 

X{xt)=Xo{xt) + -^Xi{xt) (5.2) 
47r 

with the leading contribution Xq{x) resulting from Z penguin diagrams and box-diagrams and 
Xi{xt) denoting QCD corrections to these diagrams that will be discussed below. 
In the case of charm contributions it is useful to define the parameter 

PciX) = ^ (^X^(xc) + Ix^i^c)^ , (5.3) 
with A = l^^sl being the Wolfenstein parameter (A ~ 0.225). 
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Keeping terms to first order in as, the perturbative expansion of Pd^) fias the following 
general structure 

Pc{X) = ^Pi°)(X) + P('HX) + ^PPiX) , (5.4) 
asific) 47r 

In the case of the decay Kj^ — t- tt^i^u only the top function X[xt) matters. Similarly in the 
case of S — >• Xsyi> only this function matters. 

5.2 Effective Hamiltonians for — )• and Bs^d 1^^ 

In the case of Kl — t- fi^ only subleading short distance (SD) part can be computed. The 
analysis of this part proceeds in essentially the same manner as for — )■ ti^vv. The only 
difference is introduced through the reversed lepton line in the box contribution. In particular 
there is no lepton mass dependence, since only massless neutrinos appear as virtual leptons in 
the box diagram. 

The effective hamiltonian in next-to-leading order can be written as follows: 

G OL 

HeffiKL ^ = — ^— {V:sVcdY{x,) + Vt*sVtdY{xt)) {sd)v-A{Mv-A + h.c. 

V2 27rsm Pw 

(5.5) 

The function Y(x) is given by 

Y{xt) = Yo{xt) + ^Y,{xt), (5.6) 

and 

Pc{Y) = ^ (5.7) 

has an expansion similar to Pd^) in ()5.4p . 

Only the function Y{xt) is relevant for Bg^d A*^/^^ for which the effective Hamiltonian 
reads 

n,s{Bs ^ fi+fi-) = ."2. VtlVtsY{xt)ibs)v-Aill)v-A + h.c. (5.8) 



V2 2TT 



sm^ c^w 



with s replaced by d in the case of Bd — )• 

We give here only the diagrams contributing to Yo(xt) to emphasize that only the box 
diagrams have to be calculated relative to the vu case as the direction of the internal lepton line 
differs in this case (compare with Fig. [9|). 

We are now ready to discuss NLO and NNLO QCD corrections to these decays. 
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Figure 8: One loop diagrams contributing to rare decays with charged leptons in the final state. 
5.3 NLO and NNLO QCD Calculations 

The story of NLO QCD corrections to rare K and B decays begins in the fall of 1990 when I 
started the calculation of 0{as) corrections to the flavour changing Z^'-penguin one-loop dia- 
grams that dominate semi-leptonic rare decays like K — )• vri^i/, Bg i^ — )• and B — )• Xgduu. 
The calculation involves 30 two-loop diagrams with three examples shown in Fig. [9l Most of 
them are free from infrared divergencies so that the external masses and momenta can be ne- 
glected. The infrared divergent diagrams can be regulated by non-vanishing external masses but 
a more elegant method is the dimensional regularization. 

I have started this climb by myself but after roughly twenty two-loop diagrams I stopped. I 
found this solo climb doable but I was already involved in the calculation of two-loop anomalous 
dimensions of penguin operators described previously and moreover it is always safer to have a 
partner in the climbs of that sort. Fortunately, in contrast to the ordinary climbing, in situations 
like that there is no need to return to the base camp. With good notes the climb can be continued 
whenever one decides to do it. 

In July 1991 Gerhard Buchalla returned from the military service, fit to begin his PhD 
studies and to join me in the Z*^-penguin climb. However, as a warming-up I suggested to him 
a one- loop calculation: 0{as) corrections to non-leptonic c-quark decay at ;U = 0{mc) in the 
NDR and HV schemes that in 1981 was done by Altarelli et al. in the DRED scheme. Gerhard 
found the compatibility of his results with those of O E] , published them |106| 0, and started 
from the base camp to climb the Z'^-penguin NLO QCD summit in the early summer of 1992. 
This clearly motivated me to continue my climb of 1990. Gerhard was one of my best PhD 
students ever and even if we were climbing separately, I knew that after reaching the summit, I 
would meet him there to compare my results with his. 

As a CERN fellow from 1975 to 1977 I had the opportunity to talk to one of the old masters 

" Such calculations have been refined in the context of b — > ccs |107U108] with the participation of my assistant, 
Patricia Ball. 
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Figure 9: Examples of two-loop diagrams contributing to Xi{xt)- 

of higher order calculations of g — 2, A. Peterman. He told me that in order to be sure that a 
result of a lengthy multi-loop calculation is correct, the climbing partners, should have as little 
contact with each other as possible, comparing their results only at the end of the climb. 

While in most calculations, I have done in the 1990's, I followed Peterman's advice as much 
as possible, the calculation of 0{as) corrections to Z'^-penguin diagrams with Gerhard could be 
considered as a perfect example of such an approach. Our calculations were totally independent. 
I have no idea how he got the final result and the same applies to him with respect to my 
calculation. In the fall of 1992 we compared our results diagram by diagram reaching full 
agreement on all 30 diagrams except for one term in one diagram, that turned out to be a 
misprint in my notes. We were rather confident that our result was correct. 

As a byproduct we could extract from our 0{as) corrections to the Zbs vertex the 0{as) 
corrections to the flavour conserving Z^bb penguin diagram in the large mj limit, being in fact 
the flrst group that confirmed the results of [159] done in the context of electroweak precision 
studies three months earlier. In the flrst years after the appearance of our paper |48j . it got 
most citations precisely for this additional calculation. After the discovery of — )• vr^z^P in 
1997 things of course changed. 

In the following paper [49j, after calculating 0{as) QCD corrections to AF = 1 box diagrams, 
we could flnally present 0{as) corrections to all rare K and B decays dominated by internal top 
quark exchanges: Kl — tt^vu, Bg^d — ^ /f^^^^ and B — t- ^z^P. In the case of — t- tt^vu and 
Kl — )• jU^//~ we had still to calculate NLO QCD corrections to the internal charm contributions. 
This was the subject of our third and final paper of this period [51] which required the inclusion 
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of QCD corrections to the bilocal structures as the ones given in Fig. [TOj 




(a) (b) 



Figure 10: Bilocal Structures. 

Our calculations of two first papers reduced the 0(15%) uncertainty in the branching ratios 
for Kl — )• vr'^z/P, B^.d — ^ /^^M^ and B — )• Xg^dt^i? due to the choice of fit in int{fit) present in 
the LO calculations in [U] down to ±1%. Amusingly the most important phenomenological 
contribution of the second paper is the realization that most, if not all, papers in the literature 
missed an overall factor of 2 in the branching ratio for Bg — )• fi^ fi^ . When I mentioned during 
an experimental discussion at Beauty 1995 in Oxford that the simulations for the LHCb for 
this decay should use a branching ratio by a factor of two higher, there was a real joy among 
my experimental colleagues. Afterall a branching ratio of 3.2 x 10~^ is easier to measure than 
1.6 X 10-9. 

Concerning QCD calculations the most important achievement of these three papers is the 
the calculation of Pc^\x) in ()5.4p . This allowed the reduction of the uncertainty due to the 
choice of fic in mdpc) in the -fc(^) from ±26% in the LO down to ±10%. Ten years later this 
calculation has been extended to the NNLO level, the last term in (|5.4|) . by Martin Gorbahn, 
Ulrich Nierste, Ulrich Haisch, still another multi-loop star among my PhD students, and my- 
self [31], reducing the uncertainty in question down to ±2%. This should be a very relevant 
improvement when data on — t- vr^i^P will become accurate. As in the most extensions of 
the SM, the charm contribution to — )• vr+z^z/ remains essentially unaffected by new physics, 
these results are also relevant for most extensions of the SM. 

In the case of Kl — )• fi'^ fi~ our NLO calculation reduced the /Xc uncertainty in PciX) from 
■iiAA% present in the LO down to ±22%. Ten years later this calculation has been extended to 
the NNLO level by Martin Gorbahn and Ulrich Haisch [52] reducing this uncertainty down to 
±7%. As the charm contribution is much less relevant in this decay than in — )• tt^vv and 
Kl — >■ /U^/x^ is subject to non-perturbative uncertainties, this left over uncertainty is practically 
negligible. Finally, our calculations of NLO QCD corrections to Y{xt) basically eliminated fit 
uncertainty both in Kl — )• fi~^fi^ and Bg d — )• /i^/x". 
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In September 1998, almost five years after the completion of our Z -trilogy, that also included 
QCD corrections to the relevant AF = 1 box diagrams, I spent one month in the CERN theory 
group developing a general parametrization for B — )■ ttK decays in collaboration with Robert 
Fleischer and rather decoupled from K — t- ttvu decays. Also Gerhard Buchalla was there. One 
day we received an e-mail from Mikolaj Misiak who in collaboration with Jorg Urban calculated 
0{as) corrections to the top contributions to rare decays |160j . the subject of our first two 
papers. Their result for Z^-penguin contributions obtained using the external masses of quarks 
as infrared regulator agreed with our result but the one for box diagrams disagreed with our 
calculation. The difference was phenomenologically irrelevant (1—2%) but there was a difference. 
I was already worried that I had to go again through my lengthy two-loop calculations of 1992 
but fortunately we got a second e-mail. Mikolaj and Jorg did also the calculation of box diagrams 
regulating infrared divergences dimensionally, as we did five years before, confirming our results 
in the full theory. This was a true relief. However, they found that our very simple calculation 
in the effective theory did not include an evanescent operator related to the dimensional infrared 
regulator that had to be added in the case of the box part. After including it, the small difference 
between our results disappeared. 

The issue of including the evanescent operators in this case has nothing to do with the 
cases discussed in section 3 and is rather sophisticated. It is discussed in details in [1601 1161] . 
This example shows that the method of using dimensional regularization to regulate infrared 
divergences and not distinguishing this regulator from the dimensional regulator of ultraviolet 
divergences, although very elegant and correct, does not allow a good test of the final result. In 
this respect regulating the infrared divergences by external quark masses, as was done in our 
calculations of r]2 and r]B and also by Mikolaj and Jorg in their calculation of rare decays, is 
more difficult but safer. 

Mikolaj and Jorg checked only our calculation of top quark contributions identifying the 
small difference mentioned above but somehow they were not interested in looking again at 
the internal charm contributions. This is what Gerhard and I did in |161] adding the small 
contribution of the evanescent operator to our old result. Our 1999 paper can be considered as 
a compendium of all expressions for the Wilson coefficients relevant for rare decays K — )• iri'iy, 
Bs,d ~^ f^^l^^ Slid B — )• Xs^dT^v in the SM at the NLO level. The hard work has been done in 
our first three papers that are also summarized and discussed in our review [I5j but the final 
fully correct expressions at the NLO level are given in our 1999 paper. 

The NNLO QCD calculations for the charm component in — t- ti'^vU [3T] and in Kl — 
/i [52] were of course much more involved and included several two-loop and in particular 
three loop diagrams in the bilocal operator sector. The result was the Pi '{X) in the case of 
— )• tt'^vv and the Pc iY) in the case of Kl — )• /i"*"^". I will not describe them here as 
they have been described in great detail in two papers on — )• tt^vu I was involved in and 
in the paper by Gorbahn and Haisch [52] on Kl — > fi^n^- These calculations further reduced 
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the perturbative uncertainties in these decays to the level of ±2% so that the most important 
uncertainties in the corresponding branching ratios reside in the CKM element \Vcb\ that enters 
the branching ratios as 

In Tabled] we collect references to QCD calculations for rare K and B decays. A given entry 
means that full NLO or NNLO corrections to the decay in question have been calculated in the 
quoted paper. 



Table 4: Rare K and B decays 



Decay 


NLO 


NNLO 


Kl -n^vv, B l+r, B X^uu 

— )• TT^UD 

—7- vr^/i/i 
EW to Charm in — )• ■n^uv 
EW to Top m.K ^ txvv 


|48l l49U16nU161| 

[5ll|161j 

[5111161] 

)165l 

[166] 

[Ml [57] 


m 

[52] 



5.4 Two-loop Electroweak Contributions 

The rare decays like K — t- tivv are theoretically very clean as all non-perturbative effects in- 
vestigated by a number of authors [1621 11631 1164] have been found very small and definitely 
below any experimental sensitivity in the coming ten years. For this reason it is of interest to 
investigate also two-loop electroweak contributions to rare decays. 

At first sight these contributions appear to be negligible. This however is not fully true for the 
following reason. The effective Hamiltonian for K — )• t^vv decays involves electroweak parameters 
like Gp, o: and in particular sin^ 6w that all depend on the renormalization scheme used in the 
usual electroweak precisions studies. This dependence can only be reduced by including higher 
order electroweak corrections to the leading one- loop diagrams in rare decays. This means the 
calculation of two- loop electroweak diagrams. Consider for instance sin^ Ow- The Particle Data 
Group gives two values for this parameter: (sin'^ ^w)jis ~ 0-231 and (sin^ ^VF)efr = 0.224, and 
of course there are other possibilities. As Br{K~^ — ?• vr+i/p) is inversely proportional to sin^^vK 
the two choices give two values for Br{K^ — )• tt^uu) that differ by 6%. This is clearly irrelevant 
today but I hope that the experimental data will improve in this decade to the extent that this 
difference will matter. 

A calculation of two-loop electroweak effects is clearly a very difficult affair but fortunately 
Gerhard and me could remove the ambiguity in question approximately without doing this calcu- 
lation at all [56] • We noticed that the calculations of similar effects in the context of electroweak 
precision studies contained sufficient information to find two-loop electroweak contributions to 
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K — )• vrz/P in the large mt limit without performing any loop calculations. Adding these con- 
tributions to our previous result reduced the ambiguity in question below 1%. Moreover, with 
the choice of (sin^ 6w)j^ made by us not fully accidentally in 1994, the two-loop electroweak 
corrections to K ^ nvv' can be safely neglected and similarly for other decays considered in this 
section. 

Much harder calculation has been done by Paolo Gambino, Axel Kwiatkowski and Nicolas 
Pott [58j. They calculated full two- loop electroweak contributions to — mixing, finding 
also in this case a very small effect in the MS scheme. 

Coming back to my paper with Gerhard Buchalla on two-loop electroweak corrections to 
rare K and B decays [56], we warned the readers that our large mt limit calculation of these 
corrections could miss the true value of these corrections by a factor of two. Very recently 
the younger generation of the Munich club, Joachim Brod, Martin Gorbahn and the youngest 
member of our club Emanuel Stamou, performed full two loop electroweak calculation to if — t- 
T^uv decays, basically reaching the conclusion of our large mt limit calculation but reducing 
further the theoretical uncertainty [57] . Last but not least the electroweak contributions to the 
charm part in — )• tt^vv has been calculated by the duetto Brod and Gorbahn [166j already 
in 2008. 

Finally, I would like to note that the NLO and the NNLO summits discussed in this section 
similarly to the AF = 2 summits within the SM have only been reached by the members of the 
MNLC. However, in Sections 6 and 8 we will discuss summits which have also been conquered 
by other groups in particular those led by Christoph Greub. 

6 The B Xs-f Decay: The K2 of Weak Decays 
6.1 Preliminaries 

The calculations of NLO and recently NNLO QCD corrections to i? — )■ Xg^ decay are probably 
the best known to the physics community among all QCD calculations in the field of weak decays. 
One of the reasons is the fact that the b ^ sj transition was the first penguin mediated transition 
in B physics to be discovered in 1993 in the exclusive decay channel B — )• K*j by the CLEO 
experiment. The inclusive branching ratio has been measured in 1994 by the same group. The 
other reason is the particular structure of the QCD corrections to this decay that requires a two- 
loop calculation in order to obtain the anomalous dimension matrix in the LO approximation. 
Because of this it took six years after the first QCD calculations in ordinary perturbation theory 
to obtain the correct result for the QCD corrections to i? — )• Xgj in the renormalization group 
improved perturbation theory at LO. It involved 5 groups and 16 physicists. It is not then 
surprising that the corresponding NLO calculations took nine years. They were dominated by 
the group around Christoph Greub and by the members of the MNLC although a few other 
physicists also contributed to this enterprise as I will report below. I will concentrate here on 
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the inclusive decays as they are theoretically cleaner than the exclusive ones but the effective 
Hamiltonian is of course common to inclusive and exclusive rates. A nice review of -B — )■ Xs'^ 
including exclusive decays B — t- K*^ and B ^ pj can be found in the Flavour Bible |167j and 
in the very recent review [168]. Some comments on the exclusive radiative decays will be made 
at the end of this Section and in Section 9. 

This effective Hamiltonian for 6 — )• 57 is given at the scale fii, = 0{miy) as follows 



'V2 



Ci{flb)Qi + C7.y{flb)Q7-y + C8G{fJ'b)Q. 



8G 



(6.1) 



where in view of | V*gVub/V^*gVtb |< 0.02 we have neglected the term proportional to V*gVub- 
Here Qi....Qq are the usual four-fermion operators whose explicit form is given in (j2.5p ~ (|2.7p . 
The remaining two operators, characteristic for this decay, are the magnetic-penguins defined 
in (IXTUD . 



6.2 LO Efforts 

In 1987 two groups |169[ 1170] calculated 0{as) QCD corrections to the B — )■ Xg^ rate finding a 
huge enhancement of this rate relative to the partonic result without QCD corrections. In 1987, 
when mt < Mw was still considered, this enhancement was almost by an order of magnitude but 
with the increased value of mt in the 1990's also the partonic rate increased and the dominant 
additive QCD corrections, although still very important, amount in 2011 roughly to a factor of 
2.5. 

The additive QCD corrections in question originate in the mixing of the operator Q2 with the 
magnetic photon penguin operator Qi-y that is directly responsible for the decay b — )■ 57. The 
calculation of the relevant anomalous dimensions at LO is a two-loop affair and consequently it 
took some time before the correct result has been obtained. In 1988 Grinstein et al. [171| and 
the Canadian group |172[ 1173] calculated the renormalization group improved QCD corrections 
at LO to B ^ Xg'j using the NDR and the DRED scheme, respectively. The results disagreed 
with each other. This was clearly a surprise as the LO result for the Wilson coefficients cannot 
depend on the renormalization scheme. 

In 1990, Cella et al. [174j confirmed the NDR result of Grinstein et al., and extended it to 
include the mixing of Q2 with Qsg- The fourth calculation was done by a rising Polish star, 
Mikolaj Misiak [1761 1177 ]. who in a solo climb evaluated LO QCD corrections to i? — )• Xg'j 
decay and NLO corrections to 5 — t- Xgl^l^ decay. I will discuss the last calculation in the 
next section. In these papers Mikolaj found the fourth LO result for the decay in question and 
explained why the previous NDR calculations were incomplete. The fifth calculation was done 
by Adel and Yao jl87| who ignored the evanescent operator contributions [177] , and for this sole 
reason failed to reach the correct final result. Other papers of this period contributing to this 
discussion and formulating effective Hamiltonian for 6—7-57 transitions are [178[ I179[ . 
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Mikolaj has solved all the qualitative issues but being alone he did not succeed to reach this 
LO summit. This was done in the summer of 1993 by the Rome group led by Guido Martinelli 
with the participation of a rising Italian star, Luca Silvestrini, whose PhD thesis amounted 
precisely to solving this problem \180\ I181j . This result has been soon confirmed by the Pisa 
group around Curci \IS'2\ I183j , and subsequently by Mikolaj in an erratum to his second paper 

Each of the LO calculations preceeding the final one had its own missed subtleties. One of the 
tricky points was that while in the HV scheme the one-loop matrix elements of the QCD penguin 
operators Mi = {s'y\Qi\b) with i = 3 — 6 vanish including finite parts, this does not happen in 
the NDR scheme in which the finite pieces of M5 and Mg are different from zero. Combining 
this one-loop information with the genuine two-loop calculations of the non-diagonal elements 
of the anomalous dimension matrix involving Qi with i = 1 — 6 and Qy^ and Qsg allowed to show 
that scheme independent leading order result for the Wilson coefficients of the operators Qj-y 
and Qsg could be obtained. This important solution prompted Mikolaj Misiak to extend this 
analysis to the DRED scheme |184) and him, Stefan Pokorski, Manfred Miinz and myself |185j 
to introduce effective anomalous dimension matrices that were automatically renormalization 
scheme independent in the leading order even in the B — ?■ Xgj decay. 

Thus by the summer of 1993 the correct results for the Wilson Coefficients relevant for 
B —7- Xg'j and B — t- Xgg decays have been known at LO. However, in this year an important 
observation was made by Ahmed Ali, Christoph Greub and Thomas Mannel |186| : the LO rate 
for B — )• Xgj exhibited a very large renormalization scale dependence. Changing the scale fj.^ in 
the Wilson coefficient from mf,/2 to 2mf, changed the rate of — )■ Xgj by roughly 60% making 
a detailed comparison of theory with experiment impossible. In 1993 this was not yet a problem 
as the inclusive rate was unknown experimentally at that time but the discovery of the decay 
B —7- K*^ by the CLEO collaboration in the summer of 1993 was a signal that an inclusive rate 
will be known soon as well. In any case this problem applies to i? — >• K*'y as well. 

The large Hb dependence found at LO in this decay is actually not surprising. Afterall the 
QCD effects in this decay are very large which can be traced back, at least in part, to the large 
anomalous dimensions of the dipole operators. 

In the summer of 1993 motivated by the work of Ali, Greub and Mannel I started to look at 
the steps necessary to do a complete NLO analysis of B ^ Xgj decay with the aim to reduce 
the strong dependence found by them at LO. Manfred Miinz, my very good PhD student, 
joined me in this enterprise but we were not sure that just making an outline of this particular 
NLO calculation would be sufficient for a publication. Fortunately Mikolaj Misiak, who joined 
my group in 1993, and Stefan Pokorski, who was a visitor at the MPI for Physics at that time, 
had a complementary problem. They were investigating parametric uncertainties {as,mi„mc) 
in the B ^ Xg^ decay but similarly to us were not sure that such an analysis would be sufficient 
for a publication. Once we discovered our "problems" , it was clear that joining our efforts could 
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result in a useful paper. This turned out to be an excellent decision. Our paper |185j appeared 
in November 1993, just seven months before the CLEO's discovery of — )• Xg'^ rate and became 
the standard reference in subsequent papers in which the actual climbs of the NLO B — t- Xs^ 
summit have been done. 

6.3 NLO Efforts 

6.3.1 The Basic Structure 

The complete NLO calculation for B — )• decay consists of three rather difficult steps: 
Step 1 

The calculation of 0{as) corrections to Wilson coefficients of Qj^ and Qgg operators at 
// = 0{M]/[/). The coefficients for Qi (i = 1, ..6) at this order are known from AF = 1 non- 
leptonic Hamiltonian discussed in Section 3. 

Step 2 

The calculation of the relevant 8x8 anomalous dimension matrix at 0{a1). The 6x6 sub- 
matrix involving the operators Qi [i = 1, ..6) are known from A.F = 1 non-leptonic Hamiltonian 
discussed in Section 3. 

Step 3 

Calculation of the matrix elements {s^\Qi\h) with i = 1, ..8 in perturbation theory in as- 
In what follows I will discuss these three steps one by one in the order as given above. 

6.3.2 Wilson Coefficients of Qr-^ and Q^g at /x = 0{Mw) 

The calculation of these Wilson coefficients is much harder than the calculations of Wilson 
coefficients discussed in Section 5 due to the presence of external photons and gluons: the 
external momenta cannot be set to zero. The two-loop diagrams are then calculated with non- 
vanishing external momenta of quarks, the photon and the gluon, the result is expanded to 
second order in external momenta and masses and is matched to the result of the corresponding 
effective theory calculation. The first calculation of C^^{Mw) and C%g{Mw) at 0{as) has 
been performed by Adel and Yao already in 1993 [188]. Unfortunately it was done in the 
Zimmermann's renormalization scheme that was rather unfamiliar to phenomenologists and 
consequently in 1993 this result was not noticed by many. 

Right at the beginning of 1997, Axel Kwiatkowski (one of my assistants), my PhD student 
Nicolas Pott and myself decided to calculate these coefficients in the NDR scheme, that by 
then became the standard scheme for all NLO calculations. However, in March 1997 a paper 
by Christoph Greub and Thobias Hurth |190j appeared on the Los Alamos server, in which 
they calculated 0{as) corrections to 6*7^ (Mvi/) and C^g{Myi/) in the NDR scheme and moreover 
demonstrated that their result was compatible with the one of Adel and Yao. As by that time we 
invested already two and a half months in this calculation, this was truly an unpleasant surprise. 
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Fortunately it turned out that we still could contribute something to this calculation. Greub and 
Hurth regulated the infrared divergences in the full and effective theory by using dimensional 
regularization but to our surprise performed the matching in four dimensions. As a result of this, 
infrared divergences did not cancel in the process of matching and they had to argue that they 
will be removed by including bremsstrahlung corrections. However, they did not demonstrate 
this but only showed that if this left-over divergence was removed by hand, the resulting finite 
result in the NDR scheme was compatible with the Adel and Yao result. I am mentioning 
this here to emphasize that the inclusion of gluon bremsstrahlung effects is really unnecessary 
to obtain the correct result because on general grounds the calculation of Wilson coefficients 
can be done by choosing any external states. For this reason we continued our calculation and 
demonstrated in detail that performing the full calculation in 4 — 2e dimensions, including the 
matching, directly led to the final result of Greub and Hurth without any handwaving arguments 
for the disappearance of infrared divergences. These cancellations are very clearly seen in the 
expressions presented in |19H I192j . In any case by the summer of 1997 three groups found the 
Wilson coefficients C7^{M\y) and CsgiMw) and consequently the first step of this K2 climb was 
completed. Few months later an Italian group consisting of Marco Ciuchini, Giuseppe Degrassi, 
Paolo Gambino and Gian Giudice [U7j confirmed these result, working in the off- shell operator 
basis in contrast to the on-shell basis used by us and Greub and Hurth. 

Before continuing I would like to emphasize that in spite of my critical remarks on Greub- 
Hurth calculation in question both authors played very important roles in the study of QCD 
corrections to i? — )• Xg'j and B — )• Xgl'^l' decays. In particular Christoph Greub is one of the 
great masters of these decays and his group made important contributions here both in the SM 
and beyond it. In this context his numerous analyses of these decays within 2HDM and the 
MSSM with Francesca Borzumati, the leading lady of NLO QCD calculations in weak decays, 
should be mentioned [7S1 [UB] . 

6.3.3 Anomalous Dimension Matrix 

The anomalous dimension matrix relevant for the B — )■ Xg'j decay at the NLO level consists 
of the 6x6 two-loop mixing matrix of four-fermion operators {Qi,---Q6) discussed already in 
Section 3, the two-loop 2x2 matrix describing the evolution and mixing under renormalization 
of magnetic operators and Qsg and finally the three-loop 6x2 matrix describing the mixing 
between {Qi,...Qg) and (Q7^,Qs.g)- 

In 1994 the two-loop 2x2 and three-loop 6x2 matrices were still unknown. Mikolaj Misiak 
and my PhD student Manfred Miinz decided to perform the first of these calculations in early 
1994. While certainly not an easy task, it has been achieved by developing a new technique 
of regularizing IR divergences with a common spurious mass parameter. Such a regularization 
for gluon lines had been previously thought to be prohibited because it breaks the QCD gauge 
invariance. However, the breaking turns out to be harmless for the RGE parameter calculations 
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in mass-independent regularization schemes, so long as sub divergences are treated in a careful 
manner [196] . Understanding this fact was a milestone for subsequent calculations of beta 
functions and anomalous dimensions at the three- and four-loop levels. The paper of Mikolaj 
and Manfred appeared on the Los Alamos server in September 1994 |193] . Two months later, 
Mikolaj presented the results at the Ringberg meeting, triggering objections concerning the IR 
regularization from David Broadhurst and Kostja Chetyrkin - the great masters in difficult 
multi-loop integrals. However, Kostja realized very quickly that the method is correct, and 
offered his help in continuation of the NLO project at the three-loop level. 

The three-loop calculation of the 6x2 matrix in question can certainly be regarded as the 
most spectacular achievement in the field of higher order QCD calculations in weak decays 
in the 1990's. Kostja provided soon a very efficient recursive formula for three-loop integrals 
that allowed to begin this project. With 800 three-loop Feynman diagrams this was still a very 
complicated project led by Mikolaj with a great help from Manfred. Additional difficulty was the 
treatment of 75 at the three- loop level. The calculation in the t'Hooft-Veltman scheme would be 
simply too much time-consuming and Mikolaj decided to use the NDR scheme. Unfortunately 
my technique (see Section 3) to avoid the dangerous traces with 75 , so successful at the two- loop 
level, fails at three loops. Kostja insisted that it must be possible to define the effective theory 
in such a manner that no traces with 75 occur, as in the corresponding SM diagrams. Following 
this advice, Mikolaj replaced the standard basis of four-fermion operators {Qi, ■■■Qe) by another 
(rather complicated looking) set of operators, that allowed him, Manfred and Kostja to complete 
this project without any 75 problems. The result appeared first in [197| and the details of the 
calculation have been published in |196[ [99] . 

At this point I should mention that the calculation in question, was the first three-loop cal- 
culation in the field of weak decays. It's complexity required the use of powerful PC's and has 
been fully done by using computer programs for algebraic manipulations like [110]. Such pro- 
grams have been subsequently further developed by my PhD students, Ulrich Haisch, Christoph 
Bobeth, Martin Gorbahn, and Thorsten Ewerth and by my assistant Jorg Urban, so that similar 
techniques could be used during the last decade for the calculation of two-loop contributions in 
super symmetric theories. More about it later. 

The three-loop calculation described above together with the initial conditions discussed 
in 6.3.2 and the two-loop matrix elements of the relevant current- current operators discussed 
below allowed in 1996 for the first time an almost complete (see below) NLO analysis of the 
B —7- Xsj decay [197]. Since then all analyses of this decay used the results of 1197]. It was then 
fortunate that this result has been confirmed by Paolo Gambino, and two of my PhD students 
Martin Gorbahn and Ulrich Haisch in 2003 [194j . who subsequently extended this very tedious 
calculations to other operators as discussed in section 7 below. They confirmed also the 1994 
results of Mikolaj and Manfred for two-loop mixing of the dipole operators. 
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Table 5: i? — )• Xs"f at NLO and NLLO. 7(Mixing) stands for the mixing between 4-quark 
operators and magnetic penguins. For more references to i? — t- K*{p)^ see text. 



Decay 


NLO 


NNLO 


Ci{Mw) 


|188| 119(^11911 11921 167| 


[54i 1551 1198] 


l{Q7-t,Q8G) 
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7(Mixing) 


[1971119611194] 


1195J 


Matrix Elements 


|2UUl 12011 12U21 12031 12041 1205] 


|210]-|224] 


B K*{p)-f 


[2411 12251 1306] 


[2261 1227] 



6.3.4 Matrix Elements 

The final step of the NLO program for the B — )• X^-f decay is the calculation of the matrix 
elements (s'ylQilB). The bremsstrahlung contributions to the matrix elements of Qj-y, Qsg and 
Q2 have been calculated by Ali and Greub in 1995 [200] . This calculation has been confirmed by 
my diploma student Nicolas Pott, who extended Ali and Greub calculation to penguin operators 
[201] . In these papers also one- loop virtual contributions of the matrix elements of Q^^ and Qsg 
have been calculated. 

Much harder is the calculation of the two-loop matrix elements of the four-quark operators, 
that are relevant for the NLO analysis of -B — t- XsJ- In particular, as already stressed in our 
1993 paper [185] , the /i-dependence of two-loop matrix elements should significantly cancel the 
strong //^-dependence of the LO branching ratio, pointed out by Ali, Greub and Mannel in 1993 
[I86j . 

The difficulty in this calculation is that an expansion in external momenta cannot be made 
and one has to face a two- loop calculation with full kinematics involved. On the other hand 
in the case of the matrix element of the current-current operator Q2 an expansion in powers 
of m^/m^ can be made. The first calculation of this type, using the technique of Gegenbauer 
polynomials, has been done by Christoph Greub, Tobias Hurth and Daniel Wyler already in 
1996 [2021 1203] . As anticipated in 1993, this contribution decreased the strong /xf,-dependence 
of the rate from ±30% found by Ali, Greub and Mannel down to approximately ±5%. 

In the summer of 1997, while lecturing at the Les Houches summer school on the weak 
effective Hamiltonian and higher order QCD corrections, I started thinking about repeating the 
1996 calculation of Greub, Hurth and Wyler. As many of the members of the MNLC were 
involved already in other projects, I started looking for additional collaborators. Fortunately 
among the many very good students of this summer school, there was one who was already a 
two-loop climber at that time: Andrzej Czarnecki. Already in 1997 Andrzej had on his account 
two very important two-loop calculations [132] relevant for the extraction of the CKM element 
Vcb and knew another technique that could in principle be used for the calculation of the matrix 
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elements in question: the technique of heavy mass expansions. After a short discussion we 
agreed to join our forces as soon as we complete the projects we were both involved at that 
time. Accidentaly our discussion has been documented in the form of a photo taken by Rajan 
Gupta, that can be found in the proceedings of the school just before my contribution. 

However, in the next two years we were both very busy with other projects. Among other 
things Andrzej with his well known g — 2 calculations and I with writing up my Les Houches 
lectures [16] which thanks to great generosity of Rajan Gupta amounted finally to 250 pages. 

In 2000 Andrzej Czarnecki and me met again and decided to increase our team. We were 
joined by the master Mikolaj Misiak and by experienced, already at that time, young NLO 
climber Jorg Urban. In the first climb |204j led by Andrzej Czarnecki we confirmed the result 
of Greub, Hurth and Wyler by using the method of heavy mass expansions, generalizing their 
result to a few additional terms in the expansion in m^/m^ that was found to converge rapidly. 
In the second climb |205j led by Jorg and Mikolaj we succeeded to express all the two-loop 
matrix elements in terms of four compact two- and three-fold Feynman-parameter integrals. It 
allowed us to complete the NLO B — )• Xg^ project by calculating its last ingredient: the two-loop 
matrix elements of QCD penguin operators. This last part involved few diagrams with internal 
b-quark, implying the replacement of m^/m^ by 1 and making the expansions used in [2U2\ \'2(}'S\ 
and |204j useless for the calculation of the matrix elements in question. 

This last paper |205| has been completed in the spring of 2002, just a few months before the 
celebration of Stefan Pokorski's 60th birthday at the Architects house in Kazimierz in Poland. 
Thus eight and a half years after the outline of the NLO B — J- Xg^ program [185J, Mikolaj and 
me could summarize the result of these efforts in the volume of Acta Physica Polonica dedicated 
to Stefan's 60th birthday [206]. 

Before entering the NNLO story of this decay let me mention that very recently tree-level 
contribution b — >• uus^ to i? — t- Xgj has been calculated by Kamihski, Misiak and Poradzihski 
|207j . However, for the photon energy cutoff of Eq = 1.6 GeV or higher this contribution amounts 
to 0.4% at the level of the branching ratio. 

6.4 B Xs-y at NNLO 

The story of higher-order calculations of the B — )• Xs'j rate is not finished yet. While these 
calculations decreased the //^-dependence present in the LO expressions significantly, a new 
uncertainty has been pointed out by Paolo Gambino and Mikolaj Misiak in 2001 [208j . It turns 
out that the B — )• rate suffers at the NLO from a significant, ±6%, uncertainty due to the 
choice of the charm quark mass in the two-loop matrix elements of the four quark operators, in 
particular in {s^\Q2\B). 

In the first calculations [2021 12031 12041 1205] the pole charm quark mass has been used. But 
as stressed by Mikolaj and Paolo there is no particular reason why the pole mass should be used 
instead of the MS mass, and in fact they argued that the latter choice is more appropriate in 



62 



the case at hand as charm appears only as internal particle. 

While the arguments of Mikolaj and Paolo are very plausible, it is clear that finally the 
B —7- Xg'j rate cannot depend on the choice of the charm quark mass, even if this dependence 
is significant at the NLO level. 

Thus in order to remove or reduce this uncertainty a NNLO calculation is necessary, a truly 
formidable task. It requires various calculations in three steps: 

Step 1: 

C7^{Mw) and CsgiMw) at the three- loop level and those of Ci{My/) (i = 1 — 6) at two- loop 
level, 

Step 2: 

Three-loop 6x6 and 2x2 anomalous dimension matrices of the operators {Qi, ■■Qq) and 
(Q77, Qsg) as well as the mixing between these two sets of operators at the four-loop level! 
Step 3: 

The matrix elements {s^\Qi\B) (i = 1, ..6) at the three loop level and of the corresponding 
matrix elements of the magnetic (dipole) operators {Qj^iQsg) at the two-loop level. 

In the spring of 2003, another workshop in the Ringberg castle took place. This workshop, 
organized by Andre Hoang, Gerhard Buchalla, Thomas Mannel and myself, gathered experts 
in heavy flavour physics and had a much broader spectrum of topics than the seminal 1988 
workshop at which the NLO story has been initiated. At this workshop Uli Haisch presented 
the first steps of the three-loop calculations of the anomalous dimensions of Qj {i = 1, ...6) 
relevant for B — )• Xg^f at the NNLO level, and of three-loop mixing between Qj (i = 1, ...6) and 
the semi-leptonic operators relevant for the decay B — )• Xgl^l^ at NNLO. We will discuss this 
last topic in Section 8. Moreover, Mikolaj Misiak outlined three-loop calculations of the relevant 
matrix elements and of the Wilson coefficients C7^{Mw) and Csg{Mw) |209j . 

In the following years after this 2003 workshop, considerable progress in the NNLO program 
of S — 7- Xsj has been made, and the B — )• Xg'j rate at NNLO can already be estimated. It was 
an effort of more than 17 theorists |228j and required a number of calculations over the period 
of six years by several groups. As I was not involved in this impressive project I will leave the 
description to the participants of this K2-like climb. I found in particular the summaries of 
Mikolaj Misiak [229], who led these efforts and of Uli Haisch [ 230J very informative and clear. I 
would like to thank Mikolaj for improving my insight in these calculations. 

Let me then only pay the tribute to the successful climbers of the K2 of weak decays by listing 
hopefully most important calculations which led the team of Mikolaj Misiak to this important 
victory. Here the summaries of Uli Haisch and Mikolaj Misiak quoted above were very helpful. 
Explicitly: 

• C7-y{M\Y) and CggiMw) at the three-loop level have been calculated by Misiak and Stein- 
hauser {198j and those of Ci{Mw) (i = 1 — 6) at two-loop level by Bobeth, Misiak and 
Urban 
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• The three-loop 2x2 anomalous dimension matrix of dipole operators has been calculated 
by Gorbahn, Haisch and Misiak |199j and the three loop 6x6 anomalous dimension 
matrix of {Qi, ■■Qe) operators by Gorbahn and Haisch [93]. Finally the four loop mixing 
of {Qi,Q2) operators with the dipole operators has been found by Czakon, Haisch and 
Misiak [195) . The latter one is the most impressive part of this grand NNLO project. It 
has involved the computation of more than 20000 four-loop diagrams and required a mere 
computing time of several months on around 100 CPU's. 

• A very difficult part of NNLO calculation turned out to be related to the last step of the 
program, the calculation of the matrix elements to the desired order. This calculation is 
not yet completed but a fantastic progress has been made as recently summarized very 
systematically by Mikolaj Misiak |,229j . 

As several groups took part in the latter step it is essential to refer to all the existing 
calculations. To this end I found it most convenient to shorten the most recent summary of 
Misiak [229]. The discussion of non-perturbative contributions is also included there. 

Once the Wilson coefficients have been calculated to the desired ordered, following Misiak 
the partonic decay rate is evaluated according to the formula 

8 

r(6^X,7),,>,„ =iV 5^ a{f,t)Cj{fit)G^j{Eo,fit), (6.2) 

where = | V^*Vtb|^ (G|.m^aem)/(327r^). At the Leading Order (LO), we have Gij = SiiSjj, 
while the 0{as) NLO contributions have been summarized above and in [206) . 

At the NNLO, it is sufficient to restrict the attention to i,j £ {1, 2, 7, 8} because the Wilson 
coefficients of QCD penguin operators 6*3^4^5^6 (/^f<) are very small. Treating the two similar 
operators Qi and Q2 as a single one (represented by Q2) we list the papers where six independent 
cases of the NNLO contributions to Gij have been calculated. 

First G77, G78 and G27 involve the photonic dipole operator Qj-y. While G77 was found 
already several years ago |21H 12121 12131 12141 1215) . the complete calculation of G78 has been 
finalized only very recently [216[ 1217) . Evaluation of G27 is still in progress (see below). 

The remaining three cases (G22, G28 and Gss) receive contributions from different classes of 
diagrams. Diagrams involving two-body final states are IR-convergent and are just products of 
the known NLO amplitudes. Three- and four-body final state contributions remain unknown 
at the NNLO beyond the BLM approximation [231) . The BLM calculation for them has been 
completed very recently [2191 1220) providing new results for G28 and G88 , and confirming the 
old ones [218) for G22. The overall NLO -|- (BLM-NNLO) contribution to the decay rate from 
three- and four-body final states in G22, G28 and G88 remains below 4% due to the phase- 
space suppression by the relatively high photon energy cut Eq. Thus, the unknown non-BLM 
effects here can hardly cause uncertainties that could be comparable to higher-order 0{a^) 
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uncertainties in the dominant terms (G77 and G27). Misiak concludes that the considered Gij 
are known sufficiently well. 

Thus the only contribution that is numerically relevant but yet unknown at the NNLO 
is 6*27. So far, it has been evaluated for arbitrary nic in the BLM approximation \221\ 1218) 
supplemented by quark mass effects in loops on the gluon lines [224j . Non-BLM terms have 
been calculated only in the nic 3> mi,/2 limit [222\ I223j . and then interpolated downwards in 
rric using BLM-based assumptions at nic = 0. Such a procedure introduces a non-negligible 
additional uncertainty to the calculation, which Misiak estimates to be at the level ±3% in the 
decay rate. 

As a first attempt to improve the situation, a calculation of G27 at rric = has been under- 
taken |232i 12331 1234j . Moreover a recently started calculation [235j for arbitrary rric is supposed 
to cross-check the nic = result and, at the same time, make it redundant, because no interpo- 
lation in nic will be necessary any more. Good luck! 

The references to -B — )• Xg^ NLO and NNLO calculations are collected in Table [5j In several 
of these papers also the less important decay B — )• X^gluon has been analyzed. We refer in 
particular to [ml EMI [Ml [551 ISMl [237] and more recent reviews [Ml MSI- Also NNLO 
calculation for B — )■ Xg'yy has been performed in [323j . 

These are truly impressive calculations and achievements. While I did not participate in this 
NNLO calculation I am very satisfied that 8 members of the Munich club and my collaborators 
took part in this grand project. These are Mikolaj Misiak, Paolo Gambino, Andrzej Czarnecki, 
Jorg Urban and four of my former PhD students: Christoph Bobeth, Martin Gorbahn, Ulrich 
Haisch and Thorsten Ewerth. 

6.5 B ^ Xdj 

The perturbative QCD corrections to the inclusive decay B — t- have been analyzed in 

[238t \2'69\ l24Uj and their structure is totally analogous to the case of the b ^ transition up 
to obvious changes in flavour indices. But as = Kife^Td 6 — )• (^7 is not small with respect 
to A( = VtbV^*^ and Ac = ^dy^d-i ^'^^ ^^^^ *° toke, into account the terms proportional to A„. 

6.6 Exclusive Radiative Decays 

The effective Hamiltonians for inclusive radiative decays apply of course to exclusive decays as 
well. Here the additional complications are hadronic uncertainties. As these are beyond the 
scope of this writing I just refer to selected papers [2251 UMl [MD (SMI [MSI [MH [MS] and recent 
reviews [1671 1168j . 
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7 Kl^ tt^Z+Z- at NLO 



7.1 Effective Hamiltonian 

In this and the following section we will discuss two well known decays: Ki, — )• ■k^I'^1~ and 
B — ?• Xgl^l^ . The reason for collecting these decays close to each other is related to the fact 
that the results for the first one are very helpful in reducing the work necessary to obtain the 
QCD corrections to the second decay. 

Thus the effective Hamiltonian for — )• T:^e'^e~ given in ()7.ip includes in addition to the 
operators Qi....Qq the semi-leptonic operators Q'jy and Qja defined in ()7.2p . On the other hand, 
the effective Hamiltonian for B — t- Xsl'^H.'^ given in (jS.ip can be considered as the generalization 
of the effective Hamiltonian for B — )• in (j6.ip to include the semi-leptonic operators Qgy 
and QiQA defined in (j8.2p . Evidently the operators in (j7.2p and (j8.2p behave identically under 
QCD interactions that are flavour blind. Only the renormalization scales differ. 

The effective Hamiltonian for — )• 'K^e^e~ at scales ^ < rric is then given as follows: 



^ [ziiii] + ryj(/i)] Qi + TyTA{Mw)Q7A 



(7.1) 



where the four-quark operators Qi,...Qq are familiar by now and the new operators Qjy and 
Qta are given by 

Qjv = {sd)v-A{ee)v , Qta = {sd)v~A{ee)A ■ (7.2) 

There are three contributions to this decay: CP conserving, directly CP-violating and indi- 
rectly CP-violating. In 1993 when Markus Lautenbacher, Mikolaj Misiak, Manfred Miinz and 
myself started looking at this decay it was not clear which of these components was the domi- 
nant one although there was some hope that the theoretically cleanest component, the directly 
CP-violating one, was the dominant one. The hope was partially based on the fact that at lowest 
order in electroweak interactions (single photon, single Z-boson or double W-boson exchange), 
this decay takes place only if the CP symmetry is violated |246 j . The CP conserving contribution 
to the amplitude comes from a two photon exchange which, although higher order in a, could 
in principle be sizable. 

Extensive work over the last 15 years on the non-perturbative CP conserving component 
and indirectly CP-violating one shows that within the SM the latter component dominates 
followed by the interference between the two CP-violating components. The non-perturbative 
part calculated by means of chiral perturbation method turns out to be smaller than these 
components. Most recent discussions of this decay including also i^L — ^ 7r'^/i+/x~ can be found 
in pT7l[M8l[M9l[250] . 
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7.2 The Structure of NLO Corrections 



Back to 1993. Our goal was to calculate the coefficients yyy and that are relevant for the 
directly CP-violating component at the NLO level. In order to see what was known at that time 
let us introduce yi through 

Vi = (7.3) 
with yi having the following structure according to PBE 

mv = Po + ^^^-^Zo{xt), (7.4) 
sm t'w 



VIA = ■ 2 a ^o(^t) (7-5) 

sm 

with Yq and Zq being two one-loop master functions in the SM. 

The next-to-leading QCD corrections to the coefficients above enter only Pq- The top depen- 
dent terms where known from the work of Fred Oilman and collaborators although they have 
written them in a different manner as Fred was involved in this field before PBE was introduced. 
Moreover Pq at LO was known as well. 

The structure of QCD corrections to this decay belongs to class 7 introduced in Section 2 
which is a bit special. In this language what was known in 1993 was the leading order represented 
by an l/a^ in Po and the top mass dependent part of the NLO represented by Iq and Zq. 

Our goal was to calculate the remaining, top mass independent, NLO corrections to Pq that 
are of 0(1) in the expansion. The most difficult part in this calculation is the evaluation of 
the two- loop anomalous dimension matrix involving the operators Qi....Qq and the semi-leptonic 
operator Qjy 0. The 6x6 submatrix was already known from the work of Munich and Rome 
groups and only mixing between these six operators and had to be calculated at two-loop 
level: as the later operator has no anomalous dimension only six elements of this matrix had to 
be calculated. The relevant diagrams are given in Fig. [TTJ After all the hard calculations that 
we performed before, this one turned out to be a relatively easy one and after a month of work 
the two-loop 7x7 matrix relevant for this decay and subsequently Pq at NLO was known. 

This calculation reduced a number of ambiguities present in leading order analyses |25H I252j 
and enhanced Pq by roughly 30%. The inclusion of NLO QCD effects made also a meaningful 
use of Oi-jj^ in this decay possible. Our paper appeared in February 1994 |253j . The two-loop 
mixing of all the four-quark operators with Q^y agreed with an earlier solo calculation of one 



**Q7A has no anomalous dimensions and in addition does not mix with the remaining operators. Therefore 
the only scale uncertainty in its Wilson coefficient originates in mt{fit) in Yo{xt). This uncertainty is practically 
removed through the inclusion of QCD corrections to this function that we have made in the context of NLO 
calculation for Kl — >■ l-i-^ f^^ in Section 5. In Kl — >■ vr'^e^e" these corrections appear first at the NNLO level. 
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(1) (2) (3) (4) 




(5) (8) (9) (10) 



Figure 11: Two-loop penguin diagrams contributing to 7^^^ , i = 1, . . . ,6. Square vertices stand 
for two types of penguin insertions. Possible left-right reflected diagrams are not shown. The 
numbering of the diagrams corresponds to the notation in Figs. |l]and[6l 



of us [177], once all the convention differences had been taken into account. Appendix B of our 
common paper |253j contains a detailed description of this issue. Important information on the 
NLO analysis of -B — )• Xsl~^l~ in |177j is contained there. I will comment on that in the next 
section. 



8 B ^ Xsl+l- at NLO and NNLO 
8.1 Effective Hamiltonian 

The rare decays B — )• Xgl^l" have been the subject of many theoretical studies in the framework 
of the SM and its extensions. Most recent reviews with a very good collection of references can 
be found in |T67l[T68] . 

The relevant effective Hamiltonian at scales /i = O(mfc) is given by 

n.sib ^ sfi+^i-) = n,s{b ^ 57) - ^VtlVtb [Cgvif^)Q9V + Cwa{Mw)Qioa] , (8.1) 
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Table 6: Kl n^l+l' and B X^l+l' at NLO and NNLO. For B K*{p)l+l- see text. 



Decay 


NLO 


NNLO 


— ^ 7r°e+e~ 


|253] 




B Xst-t~ 






Ci{Mw) 


|l77ll254j 


m 


7 (Mixing) 


117711254] 


[194] 


Matrix Elements 


[1771 1254) 


[2561 12571 12581 12591 12601 1261] 






[2621 12631 12641 12651 12661 1267] 


B Xdi+e- 


|l77ll254j 


[2681 1269] 


B K*fj.+fi- 


|27Ul 12711 1272| 


!24lll226] 



where we have neglected the term proportional to V*gVub and Ticsib — )• 57) is given in (j6.ip . In 
addition to the operators relevant for B — )• Xg'j, there are two new operators: 

Qgv = {sb)v-A{p'fJ')v , QlOA = isb)v-A{p'fJ')A ■ (8.2) 

8.2 B Xsl+l- at NLO 

As far as the last two operators are concerned the structure of QCD corrections is very similar 
to the case of '7T^e~^e~ and one can use our calculations for Ki — )• vr^e+e" stopping the 

RG evolution from Mw already at /x = 0{mh). In addition the matrix element of Qgv that 
also involves the mixing with Qi — Qq operators had to be evaluated at the NLO level. The 
latter one-loop calculation was possibly the main new achievement in this paper which I did only 
with Manfred Miinz, once we realized during a lunch in TUM-mensa in Garching that this paper 
could be completed rather fast [254] . At the same time, Mikolaj Misiak followed the instructions 
from Appendix B of our common paper |253] to remove a convention mismatch in his earlier 
NLO analysis of i? — )• Xsl^l~ [177] . It led him to a phenomenological formula for the so-called 
effective coefficients that parametrize the decay rate in question. The formula was published in 
an erratum to [177] . We compared our results prior to publication and found perfect agreement. 

While the work on — >■ n^e'^e' was harder, our results on 5 — t- Xgl^l' became the 
standard reference on NLO QCD corrections to these decays, dominatly because the data on 
B — )• XJ^l^ improved dramatically in the last 15 years, while the decay Kl — )• Tr^e~^e~ has not 
been observed yet. In fact |254] is one of the most cited papers in the MNLC. 

8.3 B Xsl+l- at NNLO 

Around the year 2000 various groups calculated NNLO corrections to — t- Xsl^l^ and B — t- 
K*l^l^ putting in particular emphasis on the reduction of various scale uncertainties present in 
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our NLO calculations. An important role in these efforts played the forward-backward asym- 
metry and the point in the invariant dimuon mass sq at which this asymmetry vanishes. The 
calculation of sq is theoretically rather clean for both decays, the feature pointed out by Burd- 
man |255) in the context of B ^ K*l^l~ . A meaningful discussion of sq begins first at the NLO 
level. From this point of view the NNLO calculations in question amount to NLO corrections 
to So reducing considerably the scale dependence of the LO result that we could have discussed 
already in 1994 but somehow did not. 

A compact summary of the NNLO calculations has been recently presented by Hurth and 
Nakao |168| . As I was not involved in these efforts I will not describe them here in detail but in 
order to be complete in references as much as possible I will very briefly summarize what has 
been done for the perturbative contributions to the decays in question. This time the review in 
|168j turned out to be very useful. 

First of all let us stress that for the NNLO calculations of S — )• Xsl~^l~ many parts of 
the NLO calculations of i? — )• Xg'^ can be used. This is evident when one compares the formal 
structure of QCD corrections given in (|2.48p for Class 6 relevant for B — t- Xs'^ with the structure 
of QCD corrections given in (|2.49p for Class 7 to which B — )• Xsl~^l~ belongs. Indeed, whereas 
0{as) corrections in the latter decay represent NNLO term, they represent the NLO term in 
the case of — t- Xg'^. Therefore below I will only list the calculations that are specific to 
B — )• Xgl'^l" and were not done in the context of B ^ Xg^. Again we divide the calculations 
in three steps as in previous Sections: 

Step 1: 

In [5l] 0{as) corrections to the Wilson coefficient of Qgy have been completed. More 
explicitly such corrections to the penguin function Z were still missing, while those to the 
function Y were already known from the calculations for Bg — )• /^^/^i^ described before. In 
this manner the large matching scale uncertainty of 16% at the NLO level has been practically 
eliminated. Note that the coefficient of QioA at this order, represented by the function Y was 
also known as already mentioned in a footnote few pages before. 

Step 2: 

The 0{ag) corrections to the term Pq have been calculated by Gambino, Gorbahn and Haisch 
[194| who first generalized my two- loop calculation of the mixing between [Qi — Qq) operators 
with the semi-leptonic operator Qgy^ done in collaboration with Lautenbacher, Misiak and Miinz 
(7^6'' in (|2.37p ) to next order by calculating 7^6^. At this order also three loop mixing in the 
{Qi — Qe) sector is required. It has been calculated by Gorbahn and Haisch [53] and entered 
already B Xg'j at the NNLO level. 

Step 3: 

Let me note that until this point the calculations discussed in this section were fully in the 
domain of MNLC. However, in Step 3 at the NNLO level other groups dominated the NNLO 
calculations, in particular Christoph Greub and his powerful army. In fact the four-quark 
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matrix elements including the corresponding bremsstrahlung contributions have been calculated 
for the low-q^ region in [2561 1257} 1258] . bremsstrahlung contribution for the forward-backward 
asymmetry in i? — )■ Xsl^l~ in [2591 12601 1261j . and the four-quark matrix elements in the high- 

region in [2581 1262j . The two-loop matrix element of the operator Qgv has been estimated 
using the corresponding result in the decay mode B — )• Xulu and also Pade approximation 
methods |263j : this estimate has been further improved in [266j . For QED corrections we refer 
to [263 [2M1 [2651 IMS]. The SCET methods at NNLO have been used in [267j . 

For the study of the decay B — )• Xcil^l~ we refer to j268l I269j . The analysis is very similar 
but this time CKM suppressed terms have to be kept as in the case of i? — t- Xd^y decay. 

Finally I would like to mention a paper on NNLO corrections to i? — )• Xgl'^l' in the MSSM. 
This is the work done in collaboration with my PhD students Christoph Bobeth and Thorsten 
Ewerth [81]. The main motivation for this work was the reduction of scale uncertainties in 
So which differs a bit from the one in the SM and in order to feel this difference perturbative 
uncertainties have to be under control. I should emphasize that in contrast to the rest of 
my papers described above, my contribution to this paper was minor. Indeed Christoph and 
Thorsten should be fully credited for this work. 

8.4 B ^ KI+ 1- and B ^ K*{p)l+l- 

The NLO and NNLO QCD corrections discussed for inclusive decays can of course be also 
used for corresponding exclusive decays that are easier to measure. Again as in the case of 
B — )• K*{p)^ formfactor uncertainties matter. In particular B — )• K* jj.'^ p~ decay has been 
investigated by many authors of whom we can cite only a few. Among the older papers let 
me just mention [2701 12711 1272j . In particular in 1999, Ali et al. calculated the dilepton mass 
spectrum and in the SM and various SUSY scenarios using naive factorization and QCD sum 
rules on the light cone [271j . Later it was shown by Beneke et al. |24H 1226] that B — t- K* p^ p^ 
admits a systematic theoretical description using QCD factorization in the heavy quark limit 
— )• oo. This limit is relevant for small invariant dilepton masses and reduces the number of 
independent form factors from 7 to 2. Spectator effects, neglected in naive factorization, also 
become calculable. In [273] . a calculation of i? — )• K* p^ jjL~ using soft-collinear theory (SCET) 
was presented. 

There are other aspects related to formfactors but this topic is outside the scope of this 
review and I refer only to most recent papers [2741 12751 1276^ [277] . where the formfactor issues are 
discussed in detail and various symmetries and asymmetries in the SM and various New Physics 
models have been studied. In these papers a definite progress on the calculation of formfactors 
has been achieved. A lot of information can also be found in recent reviews |1671I168] . Finally an 
important paper discussing theoretical aspects oi B ^ Kl^l~ and B — )• K*l^l^ at large (f' , in 
particular in connection with hadron-quark duality is a very recent paper by Beylich, Buchalla 
and Feldmann |278] . 
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The next pages of this review are dedicated to the NLO and NNLO QCD corrections done 
in the QCD factorization approach to two-body B decays. This means Gerhard Buchalla is 
entering the scene and I can collect the energy for the final Section. 

9 QCD Factorization for Exclusive B Decays (by G. Buchalla) 

9.1 Introduction 

The formulation of factorization theorems for exclusive hadronic S-meson decays in 1999 made 
an entire new class of processes accessible to systematic calculations of higher-order corrections 
in QCD [279\ I28U] . These processes include B decays into a pair of light mesons, the prototype 
of which is -B — )• tttt, but also rare and radiative decays, such as i? — )• K*^ or B ^ K*l~^l~. 
In the heavy-quark limit, that is up to relative corrections of order Aqqd /rrib, the problem of 
computing exclusive hadronic decay amplitudes simplifies considerably. In this limit the decay 
amplitudes can be written as hard-scattering kernels, which are process dependent but pertur- 
batively calculable, multiplied by hadronic quantities such as — t- vr form factors, meson decay 
constants and light-cone wave functions, which are nonperturbative but process-independent. 
The decomposition into calculable hard contributions and universal hadronic quantities is in 
full analogy with the factorization of short-distance and long-distance terms that is the basis 
of almost any application of QCD to high-energy processes. Correspondingly the framework is 
refered to as QCD factorization for exclusive hadronic B decays, or QCD factorization for short. 
In the present section we review the factorization formula and give an overview of the NLO 
and NNLO calculations performed for exclusive B decays. We conclude with a brief outlook on 
highlights of QCD factorization and special applications to precision flavour physics. 

9.2 Factorization formula 

The matrix element of an operator Qi in the effective weak Hamiltonian for the decay of a S 
meson into a pair of light mesons M1M2 is given by [2791 l28Uj 

{MiM2\Qi\B) = F^^^^{ml) [ duTl{u)^M2{u) + (Mi o M2) 

+ f didudvTl\i,u,v)^B{i)^MM^MAy) (9.3) 

up to power corrections of order AqcD/^ifc. F^^^^^^^^ (m^ 1) are B — t- Mi^2 form factors, where 
mi^2 denote the light-meson masses, and is the light-cone distribution amplitude for the 
quark- antiquark Fock state of meson M. Here the light-cone distribution amplitudes are under- 
stood to include the decay constant of meson M in their normalization. These quantities 
define the nonperturbative input needed for the computation of the decay amplitudes in QCD 
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Figure 12: Factorization formula. 



factorization. They are simpler than the full matrix element on the l.h.s. of ()9.3p and uni- 
versal in the sense that they appear as well in many other processes, which are different from 
B — 7- M1M2. The T^{u) and T^^ {^,u,v) are the hard-scattering functions, which are process- 
specific and depend in particular on the operator Qi. They are calculable by standard methods 
in perturbative QCD. The formula (jO.Sp exhibits the factorization of the short-distance kernels Tj 
and the long-distance hadronic quantities F^^^^ and ^m- The factorization of the latter takes, 
in general, the form of a convolution over the parton momentum fractions S,, u, v £ [0,1]. A 
graphical representation of (]9.3p is given in Fig. [TJl where index j accounts for the possibility of 
more than a single B — t- Mi form factor. The second term (~ T^^) is distinguished from the first 
(~ T^) by the participation of the S-meson spectator quark in the hard interaction, indicated 
by the spectator line entering the kernel T^^. The spectator interaction requires the exchange 
of a hard gluon. T^^ starts therefore at order q^, whereas is of order unity, schematically 



The above description relies on the fact that in the two-body decay of the B meson the final- 
state particles are necessarily very energetic, with light-like four-momenta, in the rest frame of 
the B. A meson emitted from the hard interaction, such as M2 from in Fig. [T2l is then 
described by its light-cone distribution amplitude. At leading power in Aqqd /i^^b the amplitude 
is determined by the contribution from the light-cone wave function of leading twist, which 
corresponds to the simplest, two-particle Fock state. Higher Fock states give power-suppressed 
contributions and are therefore absent in the heavy-quark limit. For example, an additional 
energetic gluon, collinear to the light-like quark and anti-quark in meson M2 will generate an 
additional, far off-shell propagator when attached to the hard process T, which results in a power 
suppression. The properties of the light-cone wave functions, which vanish at the endpoints 
(m = 0, 1), also imply the suppression of highly asymmetric configurations where one parton 
carries almost the entire meson momentum and the other parton is soft. 





(9.4) 
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To leading order in QCD, at 0{a^), the factorized matrix element in (j9.3p reduces to a 
particularly simple result. The second term T^^ is absent at this order and T\u) becomes 
a u-independent constant. Taking the matrix element of operator Qi = {uh)v-A{du)v-A for 
B —7- 7r"'"7r~ as an exampl^, the factorization formula then states that 



{T:+Ti-\{ub)v-A{du)v-A\B) = {T:+\{uh)v-A\B) (7r-|((iw)y_^|0) = iF^^'^(0)/^m| (9.5) 



This corresponds to the prescription of factorizing the matrix element of the 4-quark operator 
into a product of matrix elements of bilinear quark currents. Such an ansatz, which has a 
long history in phenomenological applications |290j . thus receives its proper justification in the 
context of QCD factorization. The approximation in (j9.5p means that the emission of the vr" 
is independent of the remaining 13 — )■ vr^ transition. The intuitive argument for this, namely 
that the energetic and highly collinear, colour-singlet ud pair forming vr" has little interaction 
with the rest of the process, has been described in |291j . The factorization theorem (j9.3p is the 
formal implementation of this idea and it allows us to compute corrections systematically. 

Factorization also works for decays of the type B — )• D^ir^ with a heavy and a light meson 
in the final state, if it is the light meson that is emitted from the hard interaction (meson M2 
in Fig. I12p . In this case spectator scattering is power suppressed and the factorization theorem 
in the heavy-quark limit takes the form 



The expression in ()9.6p had already been used in |282| to compute the order-a^ corrections to 
the ratio of the B — )• Dtt^ and B — )■ D*tt^ decay rates, prior to the systematic development of 
QCD factorization. 

The factorization theorem can be formulated using soft-collinear effective theory (SCET) 
[2921 1293j . This formalism is useful for proving factorization |294j and for disentangling the hard 
and hard-collinear scale in explicit terms. QCD factorization and SCET are theoretical concepts 
that are fully compatible with each other, but they refer to different aspects of the problem of 
-B-decay matrix elements. In some sense the relation between QCD factorization and SCET 
is similar to the relation between the heavy-quark expansion (HQE) and heavy-quark effective 
theory (HQET) in their application to inclusive B decays. QCD factorization [279tl28U"] refers to 
the separation of the matrix elements into simpler long-distance quantities and calculable hard 
interactions, where the long-distance form factors are defined in full QCD. SCET, on the other 
hand, is a general effective field theory formulation for the relevant QCD modes (hard, hard- 
collinear, collinear, soft) and allows a further separation of scales, for instance in the transition 
form factors. However, working with form factors in full QCD often seems preferable in practice. 

^^In the present section the numbering of operators Q1.2 and their coefficients is reversed with respect to the 
notation of the previous sections. 




(9.6) 
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Figure 13: Order Os corrections to the hard scattering kernels T/ (first two rows) and T/^ (last 
row). In the case of T/, the spectator quark does not participate in the hard interaction and is 
not drawn. 



9.3 NLO calculations 

Next-to-leading-logarithmic accuracy (NLO) in the decay amplitudes requires the calculation of 
the 0{as) terms in the factorization formula (|9.3p . The relevant diagrams for the kernels T^'^^ 
are shown in Fig. [131 They provide a concrete illustration of the schematic picture in Fig. [T2l As 
an example, we quote the contribution from current-current operators Qi^2 to the B — )■ vr^vr^ 
amplitude. (The penguin contributions are given in |279l 1283] .) Up to power corrections and to 
NLO precision this amplitude reads 



'"(0)Am| [aij + aiji] + ... 



where 



ai,i 



^ C2 C2 Cpa-s 



-121nii 



18 + 



du 3 



1 - 2u 



1 



■ Inu 



(9.7) 



(9.8) 



and {u = l 



aiji 



U, V = 1 — v) 

IV N m|FS^^(0) 



f'fmo /"^<f.w t'^-oA^) 

Jo ^ Jo U Jo V 



(9.9) 



Here Cp = {N'^ — 1)/{2N) and is the number of colours. 

The contribution in (|9.8p represents the T^-part of the matrix element. Coefficients Ci^2 and 
Us are evaluated at a scale /j, = 0{mh). As it must be the case, the scale and scheme dependence 
of the a^-correction cancels against the corresponding dependence of the NLO coefficients up to 
terms of order a^. The constant —18 refers to the NDR scheme as defined in [95]. Whereas form 
factor and decay constant are scheme and scale independent, the pion distribution amplitude 
^>7r has such a dependence. Since it enters only at 0{a'^) in aij this is irrelevant at NLO. 

Hard-gluon exchange between the two final-state pions leads to a perturbative rescattering 
phase and thus to an imaginary part in ()9.8p . Together with the phase from penguin loops this 
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gives the formally leading contribution to the rescattering phase in the heavy-quark limit. The 
numerical value of the imaginary part has to be taken with caution because it has to compete 
with Aqcd/"^6 power corrections, which are hard to quantify. In any case, the rescattering 
phase is predicted to be suppressed, either by Us or by Aqcd/^t-?)- 

The spectator-scattering term aijj in (|9.9p is an additional, qualitatively different contri- 
bution, which first arises at order as- The scale dependent quantities C2 and as are evaluated 
at a scale fj,h = \/A//, representing the typical virtuality of the semi-hard (or, more precisely, 
hard-collinear) gluon in this process. 

The factorized structure of the amplitude in (|9.7p is in close analogy with the amplitude 
for other weak processes, for instance B-B mixing. All of them consist of certain long-distance 
quantities multiplying calculable short-distance functions. A technical complication specific to 
()9.7p is the presence of meson distribution amplitudes, whose factorization involves an integration 
over parton momentum fractions. 

The first calculation of hadronic two-body -B-decay amplitudes complete to NLO in QCD 
was performed in |279j for the three B — )• nir channels Bd — )• 7r"^7r~, B^ — )• vr'^vr'^ and B~ — t- 
7r~'ir^. The class of heavy-light final states B^ — )• D^*^^L~, with light meson L~ = vr^, p~ , 
, , . . was analyzed in detail at NLO in |28Uj (see also [281j ). A recent discussions 
of phenomenological applications of these modes can be found in |304j . The NLO calculations 
were subsequently extended to all B — t- vrvr and B — )• ttK channels |283j and eventually to 
all decays B — t- PP and B — )■ PV, where P (V) is a light pseudoscalar (vector) meson |284) . 
Decays into fiavour-singlet mesons of the type B — t- K^*'^r]^ ^ and their special properties were 
treated in [285] . The decays B — )• VV are more complicated due to the existence of different 
helicity amplitudes for the pair of vector mesons. Only decays into light vector mesons with 
longitudinal polarization are strictly calculable in QCD factorization. Early papers on this 
subject are |3U5l I286j . Comprehensive studies have been given in [2S7\ \2SS\ 1289) . where |288j 
also considers final states with axial vector mesons. 

The methods of QCD factorization can also be employed for rare and radiative B decays. In 
this case the dominant part of the amplitude comes from bilinear quark currents, whose matrix 
elements are directly given by form factors. However, the nonleptonic Hamiltonian contributes 
to the transition as well and requires a nontrivial application of the factorization formula. The 
NLO results for the exclusive decays B — )• K*^ and B ^ pj were obtained in \2A1\ \225\ I306j . 
This calculation involves the NLO Wilson coefficients for 6 — )• 57 [197J, which depend on three- 
loop anomalous dimensions, and two-loop virtual corrections to the matrix elements of local 
operators [203^ 120^ . The work of |241j included the generalization to i? — t- K*l~^l~ at moderate 
values of the dilepton mass . 
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9.4 NNLO calculations 



Important progress has been achieved in extending perturbative calculations in QCD factoriza- 
tion for charmless two-body B decays to the NNLO, including effects of order a^. In many cases 
this level of accuracy is probably below the size of uncertainties from other sources, in particular 
from power corrections. However, the explicit knowledge of NNLO corrections is of conceptual 
interest as it extends the factorization formula to the next nontrivial level in perturbation theory. 
In addition, there are quantities for which the NNLO effects are likely to be also numerically 
important. These are cases where a contribution is absent at 0(1) and thus the 0{as) term, a 
NLO contribution in the general counting scheme, is effectively the lowest order. Examples are 
strong phases relevant for direct CP violation, hard spectator scattering, or the color-suppressed 
amplitude coefficient 02, which is accidentally small at leading and next-to- leading order and 
therefore rather sensitive to NNLO effects. 

We briefly summarize the NNLO corrections that have been computed so far. These are, first, 
the 0{a1) one- loop hard-spectator interactionj^] for current-current operators \295\ \296\ I299j 
and for penguin contributions \297\ 1298) . Second, the two- loop vertex corrections (T^) have 
been addressed for the first time in |300) . where the imaginary part is computed explicitly. The 
corresponding real part has been obtained in [ 301 ] , completing the NNLO vertex corrections for 
current-current operators. These results have been confirmed in |303j . The first phenomenolog- 
ical analysis of exclusive B decays at NNLO has been presented in |302j for the tree-dominated 
— )• TTvr, TTyO and pp decays. 

The numerical impact of NNLO effects for the coefficients oi(7r7r) and a2('7r7r), which deter- 
mine the topological tree-amplitudes in i? — >• tttt, is illustrated by the following compilation from 
[M]: 

ai(7r7r) = 1.008 + [0.022 + 0.009i]/,„^ + [0.024 + 0.026i]j^o2 
-[0.012]//,„^ - [0.014 + 0.01H]77,„2 - [0.007]p 
= 1.0191°:°^! + (0.025lH{^) z ^ (9.10) 



a2(7r7r) = 0.224 - [0.174 + 0.075i]/,„^ - [0.030 + 0.048i]j^„2 
+ [0.075]//,„^, + [0.032 + 0M9i]jj^^2 + [0.045]p 
= 0.173ir^-(0.103l°:°il)z ' (9.11) 

In both equations the first line lists the leading order result together with the vertex corrections 
(/) at various orders in a^. Similarly, the second line displays the amplitude from hard-spectator 
interactions (//). It includes a model-dependent estimate of power corrections (P) from twist-3 

^"'"The hard-spectator interactions of 0(a'i) are counted as NNLO terms since they enter the decay amphtudes 
at the NNLO level. A different convention has been used in [2951 1297] . where the same corrections are refered to 
as NLO due to the fact that hard-spectator interactions first arise at 0{as). 
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contributions to the light-cone wave function of the pion. The third hne shows the full result 
together with the total uncertainty. The detailed input can be found in [301j . 

Prom (jO.lOp we see that the calculation is well under control and the uncertainties are small. 
Note that the NLO correction is suppressed due to small Wilson coefficients. In (j9.1ip the 
cancellation of leading and next-to-leading vertex contributions is clearly visible. This implies 
the dominance of the hard spectator term and the relativly large impact of NNLO effects and 
power corrections. In addition, the first inverse moment of the i?-meson distribution amplitude 
^B, which determines hard-spectator scattering as can be seen in (j9.9p . is not well known. An 
accurate prediction of the B — )■ vr'^vr'^ branching ratio, very sensitive to 02, is therefore difficult. 
The measured number appears to be higher than theoretical estimates. The agreement is better 
for B p^p° [302] . 

We conclude with a few remarks concerning rare and radiative decays. The radiative decays 
B — )• were discussed at NNLO in [227J, but some of the contributions needed to complete this 
order are still missing. Conceptual aspects of the factorization formula for B — )• Vj at higher 
orders in ag have been treated in |244ll307t [308] . The rare decays B K*l~^l~, B -)> pl~^l~ with 
order-Os corrections |24H 1226] . mentioned in 19.31 above, have formally the structure of NNLO 
processes, requiring the dominant Wilson coefficient Cg at NNLO |194j . 

9.5 Outlook 

QCD factorization can be applied to a large number of exclusive B decays and the associated 
phenomenology is very rich. We will not go into any details here, but restrict ourselves to 
a few remarks. The literature quoted throughout the present section contains much further 
information on these topics. 

Not all observables accessible in principle to a factorization calculation are equally useful 
in practice. Accurate estimates of power corrections are still beyond our control and effects of 
typically 10 - 20% are to be expected. This makes it difficult to compute direct CP asymmetries 
because these are sensitive to the relatively small strong phases, of which only the perturbative 
(though formally leading) part is calculable. However, there are many cases where the level 
of precision attainable with QCD factorization provides the basis for accurate predictions and 
flavour-physics tests. 

A first example are suitable ratios of hadronic and semileptonic rates, e.g. 
P(i?-^^-^0) 2,2, 



dr{Bd^TT+l-v)/dq^ 



2,Ti' f^\VudV\ai + (9.12) 



=0 

which are known as factorization tests. While not directly relevant for fiavour physics, they test 
QCD factorization independently of uncertainties from form factors and Vub- At present the 
precision is still experimentally limited. 

The parameter S{p^pJ^) of mixing- induced CP violation in B^, — s- PiPi is a rather clean 
quantity. The penguin amplitude is numerically small (even smaller than for B^ — s- vr^vr^) 
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and can be computed in QCD factorization with little absolute uncertainty. CKM phases may 
be extracted to within a few degrees [2871 1289j . Additional methods to constrain the penguin 
contribution also benefit from QCD factorization. Using — t- K'^K'^ and l^-spin symmetry 
should ultimately allow a precise extraction of the CKM angle 7 with a theory error of ±1° 
[289j. 

More generally, QCD factorization can be employed to estimate the size of SU{3) breaking 
in approaches that rely on flavour symmeties to determine CKM quantities from CP violation 
in hadronic B decays [309| . 

Promising observables are exclusive rare and radiative decays such as i? — )■ K*^, B ^ pj 
or i? — 7- K*l^l^ . They are dominated by form-factor terms, similar to semileptonic modes, but 
also receive (rather moderate) hadronic contributions. To those the framework of factorization 
can be successfully applied. 

In the upcoming era of precision experiments with B mesons, QCD factorization in the 
heavy quark limit, at NLO and beyond, provides us with an important tool to control theory 
predictions at a level adequate for discoveries in flavour physics. 

10 Summary 

Our story is approaching the end and we reached the summit from which the full field of QCD 
and QED corrections to weak decays can be seen. As of February 2011 NLO and NNLO QCD 
and QED corrections to most important decays are known. We collect in Table [7| the NLO 
and NNLO summits conquered by the MNLC and the names of its members who took part 
in various expeditions. In Table [8] we collect again all calculations with references to papers 
where the names of NLO and NNLO climbers outside our club can be found. These efforts 
increased significantly the accuracy of predictions of the SM for the full field of weak decays. 
This is an important step towards the indirect searches for New Physics through flavour violating 
and CP-violating processes. Personally I do not think that the calculations of the still higher 
orders of perturbation theory are really required and we should wait for the data in order to see 
what kind of new physics will be identified directly at the LHC and indirectly through flavour 
violating processes and more generally through high precision experiments. With the technology 
developed in the last twenty years calculations of higher order QCD corrections in the extensions 
of the SM selected by Nature should be straightforward even if often tedious. Most of these 
calculations have been by now automatized as exemplifled by a recent paper |31U| . where further 
references can be found. 

On the other hand still signiflcant progress in non-perturbative calculations is required in 
order to increase the power of tests of the SM through FCNC processes, in particular in the 
case of non-leptonic K, D and B decays. While in the case of i?-decays, QCDF presented in 
the previous Section allowed to make an important progress, the case of non-leptonic -ftT-decays 
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is a different story. In this context I am very curious whether direct CP violation in Kl — )• vrvr 
decays represented by the ratio e' /e is well described by the SM or not. Even 22 years after 
the seminal Ringberg 1988 workshop it was not possible to obtain the prediction for e'/e using 
lattice simulations. Such prediction is, on the other hand possible within the large N approach 
developed by Bardeen, Gerard and myself in the middle of the 1980's. The most recent analyses 
of e'/e in this framework show that e'/e within the SM, although on the lower side, agrees 
within 2(7 with the data |31H 1312] . Unfortunately strong cancellations between QCD penguins 
and electroweak penguins do not allow clear cut conclusions on how much space is still left 
for New Physics contributions to e'/e. We should hope that in this decade the answer will be 
provided by lattice groups. In this context an interesting paper appeared very recently |313] . 
Reviews of lattice results can be found in |314l 13151 1316j . 

Presently, much faster development takes place in the case of the parameter Bk- Indeed 
during the last three years an impressive progress in calculating Bk has been achieved by 
means of unquenched lattice simulations. While until recently the values in the ballpark of 
Bk = 0.724 it 0.008 it 0.028 |113] were circulated, the most recent message from RBC and 
UKQCD cohaborations [3l7] that uses domain wah QCD reads IBk = 0.749 ± 0.027. 

Interestingly this value is very close to Bk = 0.75 obtained in the large-N limit of QCD 
|3181I319] . Including 1/N corrections Bardeen, Gerard and myself [11 1320] found some indications 
for Bk < 0.75. A very recent more precise analysis of Gerard |321] puts this result on firm 
footing. Thus afterall our large-N result for Bk presented already at the 1988 Ringberg workshop 
and briefly mentioned at the beginning of this writing has been confirmed by much more precise 
lattice simulations more than 20 years later. In this context also the paper by Bijnens and Prades 
[322] should be mentioned. See [321] for more details. Whether we were just lucky remains to 
be seen. In any case the smallness of 1/A^ corrections to the large N value Bk = 0.75 should be 
better understood. 

I would like to thank all the members of the MNLC for fantastic time we had together 
conquering the NLO and NNLO summits of weak decays over so many years. In fact after 22 
years I can confidently state that my team was much bigger than the one of colonel Hunt and 
whereas he reached only the south coal on the way to the Mount Everest summit I was lucky 
to stand on many NLO and few NNLO summits. But similarly to his case the sponsors played 
also in our case a very important role. In addition to the Physics Department of Technical 
University Munich, the thanks go primarly to the German 'Bundesministerium fiir Bildung und 
Forschung' (BMBF) and Deutsche Forschung Gemeinschaft (DFG) which were both very helpful 
in supporting us for many years. Some support came also from the MPI for Physics and the 
German-Israeli Foundation. In the last years Martin Gorbahn, Joachim Brod and myself were 
also strongly supported by the Cluster of Excellence Universe in Munich and the TUM-IAS. 

Very special thanks go to Gerhard Buchalla for writing Section 9 of this review. I also 
thank Paolo Gambino and Matthias Jamin for some help in finding references. Finally, but very 
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importantly, I would like to thank Guido Martinelli for the discussion at the seminal supper in 
1988 in the Ringberg castle and the competition which definitely kept us working for several 
years. Among the members of his team my particular thanks go to Marco Ciuchini and Luca 
Silvestrini. 

A number of my friends read VI of this NLO story pointing out misprints and making 
comments on the full text. These arc Gerhard Buchalla, Svjetlana Fajfcr, Robert Fleischer, Jean- 
Marc Gerard, Jennifer Girrbach, Mikolaj Misiak, Matthias Jamin, Cecilia Tarantino, Andreas 
Weiler and Peter Weisz. I appreciate all these comments and thank them a lot. 
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Table 7: Summary of NLO and NNLO summits conquered by the members of MNLC. 



Decay 


NLO 


NNLO 


Current-Current {Qi,Q2) 

QCDP (Q3,<54,Q5,Q6) 

EWF {Qr,Qs,Q9,Qio) 
Br{B)sL 

inclusive non-leptonic decays 
B-Decays in QCDF 
Current-Current (BSM) 
Penguins (BSM) 
Semi-Leptonic {\Vch\, iKbl) 


Buras, Weisz 

Bur as , Jamin ,Lautenbacher , Weisz 

Fleischer 

Buras, Jamin, Lautcnbachcr 
Ball, Lenz, Nicrstc, Ostermaier 
Ball, Buchalla, Jamin 
Bartsch, Buchalla 
Buras, Misiak, Urban 
Buras, Misiak, Urban 
Ball, Czarnecki, Gambino 


Gorbahn, Haisch 
Gorbahn, Haisch 

Buras, Gambino, Haisch 
Bell, Pilipp 
Czarnecki, Gambino 


m 

m, VB 

m 

Tree-Lcvcl FCNC 


Herrlich, Nierste 

Buras, Jamin, Weisz; Urban 

Herrlich, Nierste 

Buras, Girrbach 

Buchalla, Lenz, Nierste 


Brod, Gorbahn 
Brod, Gorbahn 


EW to Charm in if+ ->■ n'^uv 
EW to Top inK ^ irvu 


Buchalla, Buras ; Misiak, Urban 
Buchalla, Buras 

Buchalla, Buras 
Brod, Gorbahn 

Buchalla, Buras; Brod, Gorbahn 
Stamou 


Buras, Gorbahn, Haisch 
Nierste 


B^Xs-1 






Ci{Mw) 

7(Mixing) 
Matrix Elements 
B K*{p)-f 


Buras, Kwiatkowski, Pott; Gambino 

Misiak, Miinz; 

Chetyrkin, Misiak, Miinz 

Buras, Czarnecki, Misiak, Pott, Urban 

Bosch, Buchalla 


Bobeth, Misiak, Urban 
Gorbahn, Haisch, Misiak 
Haisch, Misiak 
Ewerth; Haisch; Misiak 


B Xs£+£-: Ci{Mw) 

7(Mixing) 

Matrix Elements 


Buras, Lautenbacher, Misiak, Miinz 
Misiak; Buras, Miinz 
Buras, Lautenbacher, Misiak, Miinz 
Misiak; Buras, Miinz 


Bobeth, Misiak, Urban 
Gambino, Gorbahn, Haisch 
Bobeth, Gambino, Haisch 
Bell, Pilipp 
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Table 8: Summary of all NLO and NNLO Calculations for Weak Decays 



Decay 


NLO 


NNLO 


Current-Current {Qi,Q2) 


[5l[90| 


m 


QCD penguins (Q3, Q4, Q5, Qq) 


[MIEaEelETlEHlEH]; jiooj 


[93] 


electroweak penguins (Q7, Qs, Qq: <3io) 


iiniii96U97ii98i. iioni 


[53] 


Br{B)sL 


|1U21I1U31 110411105] 




inclusive A5 = 1 decays 


[51161 11061 11071 1108]: [109] 




Two-Body B-Decays in QCDF 


|279]-|289] 


[295]- [303] 


Current-Current (BSM) 


[301 [29] 




Penguins (BSM) 


m 




Semi-Leptonic ( Kb j Ktfe ) ^ Decays 


[121]-[126]. [1341 1135! 


[126]-[133].[136]-[139] 


'^\ 


[113] 


jl47] 


VB 


[1421 [M] 




V3 


[14411145] 


jl46] 


Tree-Level FCNC 


[1561 1157] 




ADMs BSM 


[301 [29] 






[150111511115211154] 






[1521 11531 1154] 






[481 1491 11601 1161] 






[5111161] 


[311 [52] 


— >■ Tr~^ lip, 


[165] 




EW to Charm in — )• 7r"'"z/P 


[166] 




EW to Top in vri/z^ 


[561 [57] 




B ^ x,r. Q{Mw) 


[188]-[190]. [1911 11921 167] 


[Ml EH [198] 




[19311194] 


[1991 


7(Mixing) 


[1971119611194] 


[195] 


Matrix Elements 


[2001 12011 12021 12031 12041 1205] 


[210]-[224] 


B^K%p)j 


[2411122511306] 


[2261 1227] 




]253] 




B X,i+i-: Ci{Mw) 


[1771 1254] 


m 


7(Mixing) 


[17711254] 


[194j 


Matrix Elements 


[17711254] 


[256]-[267] 


B Xd^i' 


[17711254] 


[2681 1269] 


B K*fj.+fi- 


[2701127111272] 


[24111226] 
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